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Scenarios for electroweak baryogenesis require an understanding of the effective po- 
tential at finite temperature near a first-order electroweak phase transition. Working in 
>- . Landau gauge, we present a calculation of the dominant two-loop corrections to the ring- 

•rH , 

improved one- loop potential in the formal limit g A <^ X <^ g 2 , where A is the Higgs 
self-coupling and g is the electroweak coupling. The limit A g 2 ensures that the phase 
transition is significantly first-order, and the limit § 4 <cA allows us to use high-temperature 
expansions. We find corrections from 20 to 40% at Higgs masses relevant to the bound 
computed for baryogenesis in the Minimal Standard Model. Though our numerical results 
seem to still rule out Minimal Standard Model baryogenesis, this conclusion is not airtight 
because the loop expansion is only marginal when corrections are as big as 40%. We also 
discuss why super-daisy approximations do not correctly compute these corrections. 
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1. Introduction 

Recent scenarios fl - ^ for baryogenesis via the baryon number anomaly of the Standard 
Model have stimulated a flurry of investigation into the details of the electroweak phase 
transition in the hot, early Universe. Sakharov's generic conditions for baryogenesis @ 
require (1) baryon number violation, (2) CP violation, and (3) disequilibrium. In recent 
scenarios, the source of baryon number violation is a Standard Model effect, arising be- 
cause baryon number has an electroweak anomaly and so is violated by nonperturbative 
physics. The rate of such violation is of order exp(— E sp (T)/T) where T is the temperature, 
E sp (T) ~ ira(T)/g w is the electroweak sphaleron mass, and cr(T) is the Higgs vacuum ex- 
pectation value (VEV) at high temperature. @~0 In almost all scenarios, disequilibrium 
is provided by bubble nucleation and expansion during a first-order electroweak phase 
transition. The requirement that the first-order phase transition be sufficiently strong 
then places constraints on models of the Higgs sector. In particular, it is necessary that 
baryon-number violation be turned off after the phase transition is completed; otherwise, 
the Universe will simply relax back to equilibrium, where the net baryon number is zero. 
Since the rate of baryon-number violation is exponentially sensitive to the VEV cr(T), it 
is important to study the minimum = a(T) of the finite-temperature effective poten- 
tial V((j), T) just after the phase transition. The purpose of this paper is to study the 
extent to which two-loop effects modify the results that have been found using one-loop 
approximations to the effective potential. 

In this paper, for the sake of simplicity, we shall focus on the Minimal Standard Model 
with a single Higgs doublet. Many of our calculations should be easily extendable to more 
complicated Higgs sectors. Most scenarios for electroweak baryogenesis are based on multi- 
Higgs models because the CP violation required for baryogenesis can be incorporated by 
CP violating interactions in the Higgs sector, where it directly affects classical processes 
involving the sphaleron. It has been suggested, however, that the standard CP violation of 
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the quark sector may by itself provide sufficient CP violation, so that even the the Minimal 
Standard Model, with a single Higgs doublet, could be viable. In the Minimal Standard 
Model, and using the one-loop ring-improved effective potential (which we review later), 
Dine et aU@ have established an upper bound of roughly 30 to 40 GeV on the Higgs 
boson mass mh(0) (measured at zero temperature) if baryogenesis is to occur at the weak 
phase transition. This scenario is then excluded by the experimental bound of TOh(O) > 60 
GeV. @ (These particular constraints may be evaded by models with more than one Higgs 
boson.) 

An important question for such limits is whether they are significantly modified by 
higher-order corrections. As we shall review in the next section, the loop expansion param- 
eter in this context turns out to be formally of order A/g 2 , or equivalently mh(0)/M w (0), 
where A is the Higgs self-coupling and g is the SU(2) gauge coupling. For Dine et al.'s 
upper bound on mh(0) of 30 to 40 GeV, it is not a priori clear whether mh(0)/M w (0) is 
small — it all depends on the exact numerical coefficients in the loop expansion. Indeed, 
some of the earliest studies of corrections to the one-loop potential found large corrections, 
but these studies were subsequently shown to be flawed. @ To resolve whether the bounds 
are reliable, we present an explicit computation of the dominant two-loop corrections to 
the one-loop ring-improved effective potential. Our objective is to compute the dominant 
corrections as the Higgs mass gets large, but not so large that the loop expansion breaks 
down. Specifically, we assume 

<7 4 «A«<7 2 . (1.1) 

(As we shall review in the next section, the lower limit <C A simplifies the calculation by 
justifying a high-temperature expansion.) From the numerical results of this calculation, 
we shall determine the significance of multi-loop corrections as a function of Higgs mass; 
when we find numerical corrections of order 100% or more for a given Higgs mass, we shall 
know that the loop expansion has indeed broken down. 
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The computation of such corrections has been previously attempted in the literature 
using a super-daisy approximation. Ill Unfortunately, as we discuss in Appendix C, the 
super-daisy approximation does not correctly compute the dominant corrections to the 
ring-improved one-loop potential. 

Naively, we expect that corrections to the baryogenesis bounds should be small because 
the sphaleron is a classical solution of the effective high-temperature theory. The actions 
of classical solutions (in this case corresponding to the Boltzmann exponent E sp (T)/T 
for baryon violation) are generally inversely related to loop expansion parameters. So in 
any case where baryon violation is turned off after the phase transition, we should expect 
that the loop expansion will be under control. Our results only marginally validate this 
conclusion. 

The casual reader who is interested only in our final numerical results should turn to 
figures 25 — 28 in section 8. 

In section 2, we review the one- loop ring-improved potential and discuss how to power- 
count diagrams to find the dominant 2-loop corrections for g 4 ^ X ^ g 2 . We will be led 
to adopt the formal power-counting rule A ~ g s for the rest of the paper. In section 3, 
we warm up to the problem of 2-loop thermal calculations by first concocting a scalar 
problem that is loosely analogous. We shall review the equivalence of various prescriptions 
for resummation of ring diagrams and we shall settle on one that implements resummation 
only for the effective three-dimensional theory obtained after heavy, non-static modes have 
been integrated out. In section 4, we turn to the simplest gauge theory case — the Abelian 
Higgs model. As well as computing the 2-loop potential near the phase transition, we isolate 
which contributions are important for shifting the VEV at the phase transition from its 
one-loop value. Section 5 is devoted to the contributions of chirally-coupled fermions. In 
section 6 we discuss non-Abelian theories using SU(2) as an example, and we turn to the 
Minimal Standard Model in particular in section 7. Numerical results for the size of 2-loop 
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corrections in the Minimal Standard Model are presented in section 8. Most of the results 
for diagrams contributing to the 2-loop potential are collected in Appendix A. Appendix 
B gives derivations of the high-temperature expansions of some quantities discussed in the 
main text. Finally, Appendix C contains our criticism of the super-daisy approximation 
to the effective potential. 

Throughout this paper we shall find it convenient to work exclusively in the Euclidean 
(imaginary time) formulation of thermal field theory. We shall conventionally refer to 4- 
momenta with capital letters K and to their components with lower-case letters K = 
(ko, k). All 4-momenta are Euclidean, with discrete frequencies ko = 2nnT for bosons and 
(2n + l)nT for fermions, unless stated otherwise. 

2. Power Counting and Review 

2.1. Pure scalar theory 

The classical, zero-temperature Higgs potential is of the form 



For simplicity, let us temporarily ignore the gauge and fermion sectors and review the 
effect of finite temperature in a theory with a single, real scalar field. At high temperature 
T (T » v), there is an additional contribution to the scalar mass; the effective potential 
is approximately of the form 



(2.1) 




(2.2) 



which is the same as the zero-temperature potential except that 




(2.3) 
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The addition of the thermal mass term above is responsible for symmetry restoration at 
high temperature, EHH and this approximation to the effective potential describes a 
second-order phase transition at T 2 m 24u 2 /X. 

Diagrammatically, the thermal mass term in (2"T2") arises from the quadratically di- 



vergent loop of fig. [[]. In the temperature-dependent piece of a quadratically divergent 
loop, the UV divergence is cut off at momenta of order T, and so fig. |l] is of order AT 2 for 
temperatures large compared to the scalar mass. 

The thermal mass above is simply the dominant term in a high-temperature expansion 
of the full finite-temperature 1-loop potential. The temperature-dependent piece of the 
full 1-loop potential is simply the free energy of an ideal gas of scalar particles with mass 
m 2 ((f)) = V c / 1 / (^)), which is the mass associated with fluctuations in the scalar field around 
a background field 0:0 

d 3 k 



V^(cf>,T) = V cl ^)+T J 



In < 1 — exp 



(2tt) 3 

+ (1-loop zero-temperature result) 
In the high-temperature limit T ^> m(<fi), 



(2.4) 



V (1) ((f),T) = Vd(4>) + const. + -"-m 2 (0)T 2 - -^m 3 ((f))T + 0(m 4 ) (2.5a) 

24 127r 

= cons t. + I (-u 2 + -XT 2 ) 4> 2 - — (-v 2 + -X<j> 2 ) T + 1a</> 4 
2 V 24 J v 12tx V 2 Y J 4! Y 

+ 0(m 4 ). (2.56) 

The constants above are temperature dependent but (p independent. Such constants are 
not relevant to studying the mechanics of the phase transition, and we shall generally 
ignore them. 

We have used the classical relation m 2 ((f)) = —v 2 + \ Xcj) 2 above; however, the effective 
value Q2.3p of v 2 is quite different from its classical value at high temperature (T>T C ). It 
is therefore important to make the replacement ( |2.3|) and use instead 



m 2 M) = -v 2 + ±XT 2 + 1a</> 2 , (2.6) 



which inserted into ( 2.5a ) yields the ring-improved one- loop potentialEHU * 

3/2 

T 

(2.7) 



V&faT) = const. + \ (-^ + ±XtA f - ^ (-^ + lAT^ + lA^ ) I 



This potential sums the dominant contributions of the one-loop ring (or daisy) graphs 
shown in fig. 0, where each quadratically divergent ring has been approximated by its 
high-temperature limit AT 2 /24. The substitution of z^ ff for v 1 corresponds to a resum- 
mation of the propagator as in fig. [| We shall treat this resummation more carefully and 
systematically in section 3. 

The potential ( |2.7|) appears to describe a first-order phase transition, as shown in 
fig. |J At the temperature Tq where the quadratic term vanishes, the potential is of the 
form — 6A 3 / 2 </> 3 T + cA0 4 and has a minimum at non-zero <f> ~ yXT. Just slightly above To, 
the small quadratic term will generate a second minimum at <p = 0, indicating a first-order 
transition. The symmetry-breaking minimum at the phase transition, labeled (p c in fig. [|, 
will occur at a point where all three terms of the potential (|2.7|) are the same order of 
magnitude. One easily concludes that 

4> c ~ yfXT, ~v 2 eS (T c ) ~ A 2 T 2 . (2.8) 

In fact, the phase transition in the pure scalar theory is known to be second order. ' The 
conclusions from the previous paragraph cannot be trusted because higher-loop corrections 
to the 1-loop ring-improved potential are large at 4> c . As discussed earlier, loops with 
quadratic UV divergences are of order AT 2 and are dominated by large momenta of order 



* This potential has an imaginary part at small (f>. A physical interpretation may be found in 
ref. 

' See any textbook on critical phenomena. Concerning the order of the transition in the 
gauged case, see Ref. |22[ | for an analysis more generally applicable than what we shall review 
below. 
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T. These 0(XT 2 ) contributions are all absorbed by using ring-improved propagators. UV 
convergent loops (and the non-divergent pieces of quadratically divergent ones) are instead 
dominated by their infrared behavior. In Euclidean space, this means loop momenta are 
dominated by ko = and \k\ ~ m. The dominant ko = piece of the finite-temperature 
frequency sum TJ2k gi yes such loops a linear, rather than quadratic, dependence on 
T. TT Including a factor of l/m e ff to make a dimensionless quantity, the cost of each loop is 
therefore AT/m e ff- * * Now consider the loop expansion parameter \T/m e e at the minimum 
4> c ~ \f\T corresponding to the apparent first-order phase transition of the ring-improved 
one-loop potential (|2.7|). Eqn. (|2.8|) implies that m e ff(</> c ) ~ AT, and so the loop expansion 



parameter is O(l), verifying that the ring-improved loop expansion cannot be trusted to 
distinguish between a first and second-order phase transition in this model. 

2.2. Gauge theories 

The situation is quite different when the gauge sector is included. As we now review, 
the first-order phase transition seemingly described by the ring-improved effective potential 
is associated with a small loop expansion parameter if the Higgs mass is sufficiently small. 
The phase transition is therefore indeed first order, and the ring-improved loop expansion 
is a valid tool for studying it. 

As in the pure scalar theory, there is a thermal contribution to the Higgs mass. It is 
of the form 

-u 2 ^-u 2 s = -u 2 + ag 2 T 2 . (2.9) 

where ag 2 symbolizes a linear combination of the squared couplings in the theory. In the 
Minimal Standard Model, for instance, 

a^- (A +^1 + ^ + 3^)^, (2.10) 



' ' An early discussion of this power counting may be found in ref . [|l6f . 
TTT Loops which are logarithmically UV divergent may cost a factor of (AT/m e fr) ln(T/m e s) 
We shall ignore the logarithms when power counting. 



where g2, gi, and g y are the couplings for SU(2), U(l), and the top quark Yukawa in- 
teraction. (We shall always ignore all Yukawa couplings except for the top quark. See 
section 7 for coupling normalizations.) Diagrammatically, these contributions again arise 
from quadratically divergent loops, such as in fig. [|. 

The temperature-dependence of the full finite-temperature 1-loop potential may again 
be interpreted in terms of the free energy of an ideal gas. Now we must sum contributions 
from all the various particles in the theory with masses rrii((p) induced by the background 
Higgs field <fi: 

d 3 k 



V<»(<l> t T) = V < *(<j>) + Y;±mT J ^ln|l T exp 



(2vr) 3 

+ (1-loop zero-temperature result) 



(2.11) 



where the sum is over all particle species i, rii is the number of degrees of freedom as- 
sociated with each species, and the upper (lower) sign is for bosons (fermions). The 
high-temperature limit T » rrii((p) is© 



AVi((j), T) = const. + — m?(0) T 2 - — - m 3(0) T + 0(mf), (bosons) 

24 lZn 

1 



(2.12) 



AVi((/>, T) = const. + 4gm|(0) T 



+ 0(mf). (fermions) 

As before, we shall generally ignore the (^dependent constants indicated by "const." above. 

For gauge bosons and the top quark, ((f)) is proportional to g(f> where g is the appro- 
priate coupling. Ignoring the Higgs contribution for now, the high-temperature expansion 
of the one-loop potential is then schematically of the form 

1 / O , O Q 1 



[-v 4 + agT'W - bg^T + — A</> 



(2.13) 



where for the Standard Model 



bg*(f?T-, [6M3(^) + 3M Z 3 (0)] 



T 

12tt 



;9'i + 



ll + 9lf /2 



T 

127r" 



(2.14) 



Because of the </> 3 term, the potential (|2.13| ) describes a first-order phase-transition. 



As before, we can ring-improve the one-loop potential to sum one-loop ring diagrams 
such as fig. In Euclidean space, the small, hard loops are quadratically divergent and 
contribute a thermal mass of order g 2 T 2 to the Aq polarization of the gauge fields. This is 
the usual Debye mass, and the ring improvement is implemented by incorporating it into 
M(0) for that polarization.® For instance, the W contribution to the cubic term fl2.14T ) 
changes from 

6M w(</>)-^ to [4M*(<p) + 2Ml w ( ( p)]^-, (2.15) 



where 



Ml{<j>) = \g\^\ Mf w (0) = + ^g 2 T 2 . (2.16) 



For Euclidean frequency kg = 0, which dominates the infra-red (IR) behavior of loops, 
we shall refer to the Aq polarization as the longitudinal polarization and the other two 
polarization perpendicular to the four-momentum as the transverse polarizations. The 
ring-improvement of the terms arising from the Z are more complicated due to mixing 
with the photon,© and we leave explicit formulas for section 7. 

We now estimate the order of magnitude of parameters associated with the phase 
transition. Consider the schematic form of the potential ( |2.13| ). (The ring improvement 
will not modify the following order-of-magnitude estimates.) As in the scalar case, the 
non-zero minimum at the phase transition occurs when all three terms are roughly the 
same order of magnitude, which yields 

<Pc ~ ^T, -4t ~ ^T 2 . (2.17) 

Table 1 shows the order of magnitude of several other parameters. Of particular interest 
is the loop expansion parameter for loops involving massive gauge bosons; it is order X/g 2 . 
We shall formally assume A <C g 2 (that is, <C M w at zero temperature) so that the loop 
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A ~ 2 3 




<t> 


f-NJ 


Viz / ' V 


T 


scalar mass 


m 


rs_/ 


(g 3 /VX)T 


„3/2 T 


transverse vector mass 


M 


r*j 


(a 4 /A)T 


2^ 


vector Debye mass 


y/M'l - M 2 


rsj 


2^ 




vector loop expansion parameter 


9 2 T/M 


rs_/ 






scalar loop expansion parameter 


XT/m 


r^j 


(A/2 2 ) 3 / 2 


2 3/2 




(T c -T )/T c 


rs_/ 


2 4 /A 




barrier in V between minima 


V 




(2 12 /A 3 )T 4 


2 3 T 4 


2-loop correction to potential 


(g 2 T/M)V 


f-NJ 


(2 10 /A 2 )T 4 





Table 1. Orders of magnitude of parameters in the <fi 7^ vacuum at the phase transition. 
The far right column gives the simplified power-counting rules, which assume A ~ g 3 . 
The entry for the 2-loop correction to the potential refers to the dominant 0- dependent 
corrections. 



expansion is well-behaved. (In section 3, we shall discuss in more detail why the vector 
loop expansion parameter is the relevant one for our calculation.) Furthermore, to justify 
our ubiquitous high-temperature expansion T ^> M and m, Table 1 shows that we must 
also assume g 4 <C A. 

Assuming g 4 <C A <C g 2 , our goal is to consistently compute the leading correction to 
the ring-improved one-loop potential. It is cumbersome to compare orders of magnitude 
of various corrections when there are two coupling constants g and A. Fortunately, the 
power counting can be simplified by formally taking 

A-/, (2-18) 

which is at the geometric center of the range g 4 <C A <C g 2 under consideration. This sim- 
plification seems to always order the relative size of corrections in a way that is consistent 
over this entire range of A. We henceforth always assume A ~ g s unless otherwise stated. 
The simplified power-counting is shown in the far-right column of Table 1. Note that the 
vector loop expansion parameter is g (rather than g 2 , as it would be at zero temperature) 
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and that the Debye mass and the transverse vector mass are formally the same order of 
magnitude. 

2.3. Multiple loops and IR disasters. 

In studies of high-temperature QCD, it is well-known that the perturbative computa- 
tion of the free energy breaks down at three loops because of infrared divergences associated 
with static (ko = 0), transverse gluons. HI Though the longitudinal gluons pick up a Debye 
mass of order gT at one loop, as discussed earlier, the transverse gluons do not. These 
IR divergences are interpreted as a sign that the transverse gluons must have a mass of 
order g 2 T which is not perturbatively calculable, corresponding to a finite screening length 
for static, magnetic fields in the plasma. The uncertainty in the free energy due to the 
incalculability of this effect is order g 6 T .* 

The situation is different in a spontaneously broken theory because the gauge bosons 
corresponding to broken symmetries are not massless. In pure electroweak theory, for 
example, the transverse W and Z bosons are massive in the symmetry-breaking minimum 
at the phase transition. One can compute to any order of perturbation theory in this 
minimum without infrared singularities. (A perturbative expansion is still not useful, of 
course, unless the loop expansion parameter is small.) However there is still a limit to how 
well the phase transition can be studied, because determining the temperature of the phase 
transition requires comparing the free energy of the asymmetric (p 7^ minimum with the 
symmetric </> = one, and all the familiar problems of high-temperature QCD arise in the 
symmetric minimum. The calculation of the free energy at (j) = breaks down at three 

* We remind the reader of a simple mnemonic for this result. The only scale in the effective 
three-dimensional theory, which is the source of the IR divergences, is g 2 T. So the free energy of 
the three-dimensional theory, which will be incalculable by perturbation theory, is order (g 2 T) 3 
by dimensional analysis. The relation between the 4-dimensional and 3-dimensional free energies 
is a factor of T, giving g 6 T 4 . 
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loops and has an uncertainty of order g®T 4 ; there is therefore little point in computing the 
free energy of the asymmetric vacuum to any better accuracy. At the phase transition, the 
typical size of each of the three terms in ( |2.13| ) is # 3 T 4 and the size of n-loop corrections 
are suppressed by (g 2 T/M) n ~ 1 , giving g n+2 T 4 . So there is no point computing the free 
energy beyond three loops near the asymmetric vacuum. Being unad venturous, we shall 
limit ourselves to computing 2-loop corrections. 



3. Warming Up With a Scalar Toy Model 

We now want to proceed to compute the dominant 2-loop contributions to the po- 
tential. To introduce the basic computational method, we wish to start with as simple an 
example as possible: a pure scalar theory of a single, real field. Sadly, we saw in the last 
section that the loop expansion cannot be trusted to study the phase transition in this 
model, because the loop expansion parameter AT/m e ff(</> c ) is order 1. However there is 
no reason we cannot compute the effective potential at larger values of <p where m e ff(</>) is 
larger so that the loop expansion parameter is smaller. This region of the potential has 
nothing to do with the phase transition in this model, but it provides a simple example 
of the computations and approximations that we later implement in gauged models. We 
shall take a classical Higgs potential of the form 

W) = -^V + ^V (3-1) 
and shall study the potential for (p of order 

(3.2) 

at temperatures T near the critical temperature T c . This is supposed to be analogous to the 
A ~ g 3 gauge theory discussed in the last section insofar as (1) there is a zero-temperature 
loop expansion parameter called g 2 , and (2) we examine the potential at <p ~ T (see Table 
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SCALAR TOY MODEL 


by assumption: 







T 


effective mass at <j> 


m{4>) 


rs_/ 




Debye mass 




rsj 


gT 


loop expansion parameter 


g 2 T/m{4>) 


rsj 


9 


2-loop correction to potential 






g 4 T 4 



Table 2. Orders of magnitude of parameters in the scalar toy model for ~ T and T ~ T c . 
The entry for the 2-loop correction to the potential refers to the dominant 0- dependent 
corrections. 

1). The effective scalar mass is then of order gT and the loop expansion parameter is order 
g 2 T/m ~ g, both analogous to the vector mass and vector loop expansion parameter in 
Table 1. We have named the scalar self-coupling g 2 in the toy model both to emphasize this 
analogy and to avoid confusion with our A ~ g 3 convention used in the gauged models. We 
hope that the reader's benefit from seeing the techniques first presented in a simple model 
will outweigh any confusion inherent in the temporary change of conventions. Orders of 
magnitude are summarized in Table 2. 

Let us apply our earlier review of power counting to estimate the two-loop contri- 
bution to the potential. Ignoring quadratically divergent loops, which will be absorbed 
by ring-improvement, the dominant two-loop contribution should be order g 2 T 2 m 2 , where 
there is an explicit power of g 2 '', an explicit power of T for each loop, and then m 2 by 
dimensional analysis. Since m ~ gT, the 2-loop contribution is equivalently order g 4 T 4 . 
When computing the 2-loop potential, we shall always drop any contributions smaller than 
0(g 4 T 4 ). In the region of <fi under study, the contributions dropped are of the same or- 
der as three- and higher-loop contributions, and so there is no point retaining them in a 
two-loop calculation. 
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3.1. The 1-loop result without resummation 

We shall leave the analysis of resumming rings for sections 3.3 and 3.4. For now, 
we focus on the details of the unimproved loop expansion in this model and start by 
carefully examining one- loop results. We need to keep terms to order g T , which requires 
one higher order in the high-temperature expansion ( |2.5q| ) discussed in the introduction. 
At this order the one-loop potential receives UV infinite contributions from usual zero- 
temperature divergences, and so we must confront regularization and counter-terms. We 
shall regularize all our calculations using dimensional regularization in 4 — 2e dimensions. 
Computing counter-terms to one-loop order, one finds that the bare potential expressed in 
terms of the renormalized coupling constants is 



Vbare(0) — — « ^bare^bare + ^jfi'bare^bare 

= - 1 -Z 1 v*<? + ^Z 2 g^\ 



Zi = l + ^g 2 - + 0(g 4 ), (3.4a) 
Z 2 = l + -^g 2 - + 0(g 4 ). (3.46) 

\x is the arbitrary renormalization scale. We shall use MS regularization, which, we remind 
the reader, corresponds to performing minimal subtraction (MS) and then changing scales 
to p, by the substitution 

- e 7E\V2 



(3.3) 



^ = Hw ' (3,5) 

where 7 E is Euler's constant. 

The next term in the high-temperature expansion of the one-loop potential is well- 
known (but we shall review the derivation in section 3.4): @ 

^V^^) = ^V hare (<P) + J[m(<j>)], (3.6) 
14 



J(m) = HK 2 + 



m 2 ) 



e 



(3.7) 



where 



. , 1 2^2 1 1 4 

const. H m 1 m 1 

24 12tt 64tt 2 

+ 0(m 6 /T 2 )+0(e). 



and the integral-summation sign above is short-hand for the Euclidean integration 

d S-2e k 

fc 

The sum is over ko = l-nnT for all integers n. 

The divergences in J(m) cancel against the counter-terms in Vbare (</>)• The reader 
should note that our definition of J(m) includes both the temperature-dependent and 
zero-temperature contributions.* Note also that MS regularization does not get rid of all 
factors of 7 E and ln(47r) as it would at zero temperature. The subscript ft above refers to 
the absence of resummation. 

3.2. The 2-loop result without resummation 

The 2-loop diagrams are shown in fig. [7|, where the heavy dots represent 1-loop, 
zero-temperature counter-terms. The propagators use the mass ( |3.8|) appropriate in the 
presence of the background field </>, and we shall usually refer to m 2 {<f) simply as m 2 . The 
crosses in the diagrams represent explicit factors of the background field <fi at vertices. 

Fig. 0a is the most straightforward of the 2-loop diagrams, giving a contribution to 
the effective potential of 

f'V™ = \g 2 [I{m)} 2 (3.10) 



* One often sees only the temperature-dependent piece in the literature, which has an 
m 4 ln(m 2 /T 2 ) term, which is not analytic in m 2 . It is important to realize that in the full 
result, which includes the zero-temperature contribution, the only term not analytic in m 2 is the 
m 3 T term. 
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where 

I{m) is related to J(m) of ( |3.6|) by I(m) = m _1 dJ/dm, and so its high-temperature 
expansion is 



I(m) = — T 2 mT -m 2 



1 i / A* 



e 



r.2 



+ 0(m 4 /T 2 ) 

(3.12) 



where we have now explicitly shown the leading term of order e. Though it will appear in 
results for individual diagrams, it turns out that 

' £ = ln (ft + 2 7e-21u2-2^ (3.13) 

is an unimportant constant because it will cancel in our final result. (A sketch of the 
derivation of t e may be found in section 3.4. £ is the Riemann zeta function.) Using the 
expansion ( |3.12| ) in the contribution fl3.10|) to the potential, and keeping terms only up to 
0(s 4 T 4 ), 

= COnSt " - ~ W 9VT2 K + ^ + ln (ft " " 6 

+ 0(g 5 T 4 ). 

(3.14) 

(When we later address the resummation of rings, we shall find that the first term, which 
is order g 3 T 4 , is absorbed by the resummation.) 

Fig. 0b is more interesting because it is logarithmically divergent in the three- 
dimensional theory (that is, when all loop frequencies kg are set to zero). As a result, 
it generates a logarithmic dependence on the mass m(<p), unlike the one- loop contribution 
or fig. [7]a. As we shall see when we later return to gauge theories, it is such logarithmic 
terms which are solely responsible, at the order under consideration, for modifying the 
VEV at the phase transition from its one-loop value. 
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Turning to specifics, the contribution of fig. [7|b is 



(3.15) 



where 



H(m) = 11 



4e 



H P %(P 2 + m 2 ) (Q 2 + m 2 )[(P + Q) 2 + m 2 } ' 
The high-temperature limit of H has been evaluated by Parwani,0 who finds 

"< m > = ^ r2 6 + *■ + in (ft + 2 in (ft + 2 - c " 

+ 0(m 2 ) + 0(eT 2 ) 

where 



(3.16) 



(3.17) 



(3.18) 



is a numerical constant which can be expressed in terms of double definite integrals of 
elementary functions. (See Ref. [|4j for details.) Substitution into ( |3.15|) gives 
1 



£ 48(4tt) 2 



g H 4> T 



\-^~ x At 2 



2X / T 2 

- 2 In 



2 + c H 



+ 0(/T 4 ). 

(3.19) 



Note the promised presence of In [m (</>)]. 

The two counter-term graphs figs. [7|c and 0d give 



(3.20) 



(3.21) 



Combing the results of this subsection with (|3.6|) yields the unimproved two-loop potential: 



.2e T /(2) 



1 



„ e yy = const. — 
r i 48 • 4tt 



g 2 mT 3 



+ 2 



i{(_^ (T) + ^ (T)T2 )__l_ CB! ,V 



1 



9 4 T 2 



24(4tt) 2 
1 



2 In 



r 2 



1 c B c H 



-*-m 3 T+ — 
12tt 4! 



+ 0(g 5 T 4 ) 



(3.22) 
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where 



m 2 = -z/ 2 + ^V, (3.23) 
2 3 4 / /2 2 



9 iT) = 9 32> ln (^ ,+ --'' ^ 



z, 2 (T) =z/ 2 



1 „2, f£ 



1 o9 In I —77 + 

32yr 2y It 2 ' 



-,2 



(3.25) 



As it should be, the result is invariant under the renormalization group to the order 
we have computed. We have chosen to write the answer in terms of the running couplings* 
g 2 (T) and u 2 (T). (At this order we need not worry about the anomalous scaling of 0.) The 
physical scales in this problem are m and T, and so, when evaluating ( 3.22]) in practice, we 



should choose the renormalization scale jl to avoid producing large logarithmic enhance- 
ments, ln(/i/T) or ln(/i/m), of yet higher-order corrections. Fortunately, since m ~ gT, 
m and T are not drastically different scales. The difference between g 2 (T) and g 2 (m) is 
of order g 4 lrig and is indeed small, and so p, ~ T and p, ~ m are both adequate choices. 
We chose to write ( 3.22|) in terms of g 2 (T) instead of g 2 (m), but this choice is merely 



convention and lacks physical significance. 



3.3. Resummation: method I 



How do I resum thee? Let me count the ways. 

We now want to implement the resummation of the dominant parts of ring diagrams 
by replacing m 2 ((f)) in our propagators by ?n 2 ff (</>) as defined in (|2.6| ) to include the thermal 
contribution to the mass. Parwanil0 (who has computed the sub-leading correction to 
the scalar mass at high temperatures) uses one method for doing this systematically and 
consistently. Rewrite the bare potential (|3.3|) as 



V haie (<f>) = \zi {-v 2 + ^9 2 tA <j> 2 + ^ 2e Z 2 g 2 <j> 4 - -^Z^T 2 ^ 2 , (3.26) 



* These couplings are run with the usual zero-temperature renormalization group and do not 
represent the use of any sort of temperature-dependent renormalization group equations. 
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where the thermal mass term has simply been added in and subtracted out again so that 
nothing is changed. But now interpret the first term as part of the unperturbed Lagrangian 
Cq and treat the last term as a perturbation. (We shall loosely refer to the last term as the 
thermal "counter-term," but we have not in fact changed the renormalization prescription 
from the usual zero-temperature one.) Nothing has changed if all orders of perturbation 
theory are summed. However, order by order, the g 2 T 2 pieces of quadratically-divergent 
sub-loops now cancel against diagrams involving the thermal counter-term; the new per- 
turbative expansion is controlled by the convergent- loop expansion parameter g 2 T/m. 
In the case at hand, this resummation replaces m(4>) by 



m eff (</>) = -v 2 + ^ W + l~ A 9 2 T 2 



(3.27) 



in our previous calculations of contributions to the two-loop potential^ and introduces 
the new diagram of fig. |8|. The box represents the thermal counter-term — the interaction 
generated by the last term in ( |3.26|) . Fig. |8| gives 



2e V^ = -±g 2 T 2 I(m eS ) 



48* 
= const. + 



1 2 rr3 i 1 2 2 rri2 

:9 mesT + 48(4^ meffT 



2c F 



(3.28) 



48-4tt' 
+ 0(g 5 T 4 ). 

Replacing m — > m e ff in the original 2-loop contributions ( |3.14j) , (|3.19| )- (|3.21| ) and combin- 
ing with ( |3.28| ) above now gives 



y (2) = const. + 



1 



-v\T) + ±g 2 (T)T 2 



1 



1 



4rp2 



24(4tt) 
12^ T+ 4T 



2 In 



(4tt 
T 2 



9 2 (T) + 



(47T) 



;C B g 



2 2 
2 C b9 V 



(.4 , 5^4 



(3.29) 



+ 0(g b T" 



' This substitution may be made directly in the results of sections 3T and 3^ except for the 
final formula (3.22) because the expansion ( p|) of m 2 (<f>) was used when putting the result in its 
final form. 
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Compare this to the unimproved result ( 3.22 ). As promised earlier, the 0(g 2 mT 3 ) term 



of the figure-eight diagram ( |3.14| ) has disappeared, canceled by the thermal counter-term 
diagram (|3.28 ). The original 0(g 2 mT 3 ) contribution arose from taking the sub-leading 



0(mT) contribution of one loop integral times the leading 0(T 2 ) piece of the other, and 
the latter is precisely what resummation is intended to absorb. 

3.4- Resummation: method II 

The implementation of resummation above is a little less natural for gauged theories. 
The scalar ring-diagram in fig. [I] is independent of momentum, whereas the quadratically- 
divergent diagrams of fig. [(] for gauged theories are not. Moreover, the polarization de- 
pendence of these g 2 T 2 contributions to the vector self-energy also depends on 
momentum!! and is simple only in limits such as K —> 0, where only the contribution to 
Iloo(O) is non-zero.* Having computed the leading contribution to the self-energy H(K), 
should one resum the vector propagator G (K) as l/[G ( 7 1 (X)+n(K)] or 1/[Gq 1 (K)+U(0)] 
or something else? The answer is that it doesn't matter. Resummation only affects pertur- 
bative expansions when IT cannot be treated as a perturbation to the inverse propagator 
Gq 1 . This happens only when K 2 <C T 2 , in which case Ii(K) »s n(0). 

To see more explicitly that resummation is irrelevant when IT <C Gq 1 , consider K 2 ~ 
T 2 . Then U(K) ~ g 2 T 2 and Gq X (K) ~ T 2 . The resummed propagator for K 2 ~ T 2 
can then be expanded perturbatively in powers of ITGo ~ g 2 ■ For example, the resummed 
propagator of the previous section may be expanded as** 



K 2 + m 2 + ag 2 T 2 K 2 + m 2 (K 2 + m 2 ) 2 (K 2 + m 2 ) 3 



* We emphasize that we are working in Euclidean space, and so by 11(0) we always mean 
II(A;o = 0, k — > 0) since ko is discrete. If one analytically continues ko to real-time and takes 
instead the limit H(ko — * 0, k = 0), the limit is completely different, giving the plasma mass for 
propagating waves A rather than the Debye mass for static electric potentials Ao. 

** Since m 2 ~ g 2 T 2 , we could expand the powers of m as well, but we have not bothered to 
do so. 

20 



The n-th term is order g 2n . Every term, except the first, gives contributions (inside any 
diagram) that cancel order-by- order against insertions of the thermal counter-term, as in 
fig. H So resummation for K 2 ~ T 2 has no effect on the perturbative expansion of the 
final result. (We shall see this even more explicitly in a moment by recomputing the 2-loop 
potential in our scalar toy model using a different resummation prescription.) 

The circumstance where resummation is relevant is the infrared region K 2 <(gT) 2 
where IT ~ Gq . Then all terms in ( |3.30| ) are the same order. Any resummation prescrip- 
tion which approximates IT (0) in the infrared region will do. In general, any resummation 



1/[Gq 1 {K)+V{K)] will work if (1) V{K) « 11(0) when K 2 <(gT) 2 and (2) V(K)<g 2 T 2 



for general K. 

Since K 2 <(gT) 2 implies ko = in Euclidean space, the resummation prescription for 
fco 7^ modes is irrelevant (as long as condition (2) is met). Since resummation is an infra- 
red phenomena, we find it convenient to adopt the following point of view, in the spirit of 
decoupling and the renormalization group. When computing the effective potential, first 
integrate out all of the heavy ko ^ modes (ignoring the issue of resummation) to obtain 
an effective, three-dimensional theory of the ko = modes. This will generate the thermal 
mass terms as well as other interactions induced by the heavy modes, which we compute 
to the desired order in perturbation theory. Only then do we finally integrate the ko = 
modes after deciding on a sensible partition of the effective three-dimensional Lagrangian 
into an unperturbed piece Cq, containing the thermal mass terms, and a perturbative piece. 
In the language of the previous paragraphs, this corresponds to choosing a resummed 
propagator 1/(G ~ 1 +V) with 



where n(0) is the dominant, 0(g 2 T 2 ) term of 11(0). (One could alternatively replace n(0) 
by the full 11(0) calculated to some order in perturbation theory, but perturbative changes 
to V have no effect on the perturbation expansion.) 




ko = , 
k ^0 , 



(3.31) 
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Though resumming only the ko = modes sounds a little more cumbersome than 
resumming all the modes, we find it algebraically simpler in the gauge theory case because 
l/[G'o 1 (X)+n(X)] and l/[Go 1 (K)+n(0)] turn out to have more complicated polarization 
dependence for ko 7^ than for ko = 0. We also find it conceptually simpler to apply 
resummation only to the modes which require it. 

As a paradigm for splitting calculations into heavy ko 7^ modes and light ko = 
modes, and to emphasize the fact that masses (or self-energies) may be treated perturba- 
tively for heavy modes but not for light ones, we shall briefly sketch the derivation of the 
high-temperature expansion ( |3.12| ) of I(m) from its definition (|3.11| ). 

2 / d 3 ~ 2e fc 1 2 ^ f d 3 " 2£ A; 1 

I(m)-n TJ ( 27r)3 _ 2e k2 + m2 + H 2^J (27r) 3-2e (2-rmT) 2 + k 2 + m 2 

= u 2e T / — I 2a 2e T W / — fc ^ J w 

^ 7 (2n) 3 - 2 ' k 2 + m 2 * ^o^J ( 27r ) 3_2e [(27rnT) 2 + k 2 } l+1 

rH + £ )^4VV mT 



(4tt)I V m2 

^ r(i + l) ( } UttJ C(2/ 1 + 2e) 



— T 2 mT -m 2 

12 4tt 16tt 2 



e 

00 / 9 \ 1 



1=0 



(3.32) 

where £ is the Riemann zeta function. In the first line we split the calculation into ko = 
and fco 7^ modes; in the second line we expand the heavy modes in powers of m 2 ; in 
the third we first evaluate the momentum integrals and then sum over n; and in the last 
line we take the e — > limit. (This derivation may be used to obtain the constant i e of 
( p,13|) .) The point of this example lies not in the mathematical details but in the fact that, 
in the sum over heavy modes, we are able to treat m 2 as a perturbation; the expansion 
of the integrand in m 2 for these modes is equivalent to the high-temperature expansion 
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in m 2 /T 2 . The ko = term, on the other hand, produces the only term not analytic in 
m 2 — the mT term at 0(e ). (One may now also obtain the high-temperature expansion 
( p77| ) of J{m) by using I{m) = m~ 1 dJ/dm.) 

We now demonstrate how the new method of resummation works in our scalar toy 
model. For the scalar case, the new method will seem more convoluted than the previous 
one, but we find it simpler to use when we compute in gauge theories. 

The resummation of the setting-sun diagram of fig. [7|b is shown in fig. n| The ko 7^ 



lines are marked "heavy" and the ko = lines marked "zero" . The double lines represent 
resummed propagators. This expression may be simplified by realizing that the second 
term is order g 5 T and hence ignorable; this term is of the form (leaving out e's to avoid 
clutter): 

1 4-2 T f d 3 P 1 ^ f d 3 k 1 
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4* r J (2tt) 3 (p 2 + m 2 eS ) (2tt) 3 (K 2 + m 2 )[(P + K) 2 + m 2 ] 



(3.33) 



+ 0( P yT\m 4 /T 4 ) + ---] 
= ^9 4 (p 2 T [0(m ln(T//2)) + 0(m s /T 2 ) + 0(m 5 /T s ) + ■ ■ ■} 
= 0(g 5 T 4 ln(T//2)) + 0(g 7 T 4 ) + 0(g 9 T 4 ) + ■■■. 
Note that we have implicitly used dimensional regularization to make sense of UV divergent 



integrals such as J d 3 pp n /(p 2 + m 2 ) and estimate their order as 



in 



n+l * 



We may compute heavy-loop diagrams, like the first term in fig. by relating them 



to our results from un- resummed perturbation theory. Fig. 11 shows such a relation. None 



* The reader may be concerned by this. It is possible, but cumbersome, to reorganize the 3- 
dimensional (ko = 0) integrals that we shall encounter in implementing resummation at 2-loops so 
that they are UV convergent. For instance, the difference between the resummed and unresummed 
results takes the form J d 3 pp n [(p 2 +?riefj)~ 1 — {p 2 + m2 ) which is more UV convergent. Also, 
we are only interested in the <j> dependence of the potential, and so may instead evaluate d/d(f> of 
the potential, which makes the UV behavior even more convergent. With enough care, one may 
draw the same conclusions about which terms are negligible from UV convergent integrals. 
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of the lines are resummed, and unlabeled lines represent the sum over all values of Jzq. 
The first term is the previous, unresummed result of (|3.19| ), the second term is 0(g 5 T ) 
as discussed above, and the last term is a 2-loop diagram in three dimensions. Combining 



with fig. 10, resummation of the setting-sun diagram is therefore implemented by fig. 12 



We now need the difference of two three-dimensional graphs, the last two terms of fig. 12 



which differ only by the mass used in the propagators. One may either compute the three- 
dimensional integrals directly or note that the the mass dependence in the four-dimensional 
result ( |3.19| ) can only come from the ko = term at this order. So the three-dimensional 
graphs are each 

1 . _ _ T / T 2 M 

(3.34) 



48(4tt) 



const. + 2 In ( — ^ 



(I? 

\rri- 



with m replaced by m e fi in the third term of fig. [12]. The total effect of the sum of terms 
in fig. [12| is therefore to simply replace m by m e g in the unresummed, four-dimensional 
result fling) : 



2e 



48(4tt 



4 2ei2n-i2 

~29 V T 



1 



In 



jp2 



+ 2 In 



m 



eff 



+ 2 



+ 0(/T 4 ). 

(3.35) 

The result is easily 



The resummation of the figure-eight diagram is shown in fig 
evaluated by referring to the review ( |3.32|) of the expansion of I(m). The three-dimensional 
piece of I(m) is the mT term, 



and the contribution due to heavy modes is everything else, 



(3.36) 



I h (m) =^T 2 (l + a £ 



1 



m 



64tt 2 

+ 0(m 4 /T 2 )+0(em 2 ). 
The resummation of the figure-eight is then 

-<fn, rti T 3 + 
1 



1 



+ In 



Ji2 



2c F 



(3.37) 



= const. - -j^g^sT 3 + -±— g 2 m 2 eS T 2 



48 -471' 
1 



48(4tt) 



; g 2 m 2 T 2 



e +6 £ + ln 2 



- 2c F 



(3.38) 



+ 0(g 5 T^ 
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Note that this is not simply the substitution of m — > m e fj into the unresummed result. 

Resummation of the counter-term diagrams of fig. [7]c and fig. |7|d do not change the 
unresummed results (|3.20|) and ( |3.21] ) through 0(g 4: T 4: ). The thermal-counter term of fig. |8| 
now only applies to the ko = modes, giving 



1 



48-4tt 



g 2 m eS T £ 



(3.39) 



For the one-loop result of ( |3.6|) and ( |3.7| ), resummation affects only the m 3 T term arising 
from the ko = modes: 



n 2e v^((/>) = ^v haTe ((P) + JrK0)], 



1 



J R (m) = const. + j^i 2 T 2 - —■„, 



1 



mt„T 



1 



64tt^ 



-m 



2 C p 



(3.40) 
(3.41) 



+ 0(m 6 /T 2 ). 

The sum of all these contributions reproduces the previous resummed, two-loop result 
(j3.29|) , giving an explicit example that the exact details of the resummation prescription 
are unimportant. 



4. The Abelian Higgs Model 

We now turn to the simplest example of a spontaneously-broken gauge theory: the 
Abelian Higgs model, defined by the Lagrangian 

£ = —-^F 2 + \D&\ 2 — "l/(|$| 2 ), F(|$| 2 ) = -z/ 2 |$| 2 + ^A|$| 4 , (4.1) 

where $ is a complex field and -D^$ = (<9 M — ieA M )$. We shall typically express the 
potential in terms of <p = v2$ so that it takes the canonical form (|2.1] ). We now return 
to the original power-counting rules of Table 1 and assume A ~ e 3 . 

Now consider the masses of particles in a background field <p. In most gauges, such a 
background field induces mixing between the scalar and the unphysical (/c M ) polarization of 
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the vector, as shown in fig. TA. When studying effective potentials V(<f>), mixing arises even 
in (consistently defined) gauges because there is no single, ^independent gauge-fixing 
condition that will eliminate mixing for all values of the background field </>.[0'il How- 
ever, there exists a gauge choice — Landau gauge — where the mixing effectively vanishes, 
because in this gauge the unphysical vector polarization does not propagate. In order to 
avoid having to diagonalize propagators, which increases the number and complexity of 
Feynman diagrams, we shall restrict ourselves to Landau gauge. Our gauge-fixing condition 
is therefore 

C e , f . = -^(d.A) 2 -fjd\ ^^0, (4.2) 

where the ghost r\ is completely decoupled in this Abelian model. The Euclidean vector 
propagator in Landau gauge is 

G ^(K) = K 2 +M 2 ( 4 - 3 ) 

In the best of worlds it would be nice to compute results in a variety of gauges and check 
that physical quantities are indeed gauge invariant, but we have not had the perseverance 
to do so. 

In Landau gauge, the masses of fluctuations in the background <fi are classically 

M 2 ((p) = e 2 (p 2 (vector) 

m \{4>) = -v 2 + \^4> 2 (physical Higgs) (4.4) 
m|(0) = — v 2 + |A(^> 2 (unphysical Goldstone boson) 

We can now easily construct the (un-resummed) one-loop potential: 

l? e V^ ] = At 2e T4 are (0) + (3 - 2e)J[M(0)] + J[mi (</>)] + J[m 2 (0)] + const. (4.5) 

(where the extra constant arises from the decoupled ghost contribution). As before, the 
singularities in the expansion (|3.7|) of J(m) cancel against those in the bare potential Vb are . 
The explicit counter-terms are given in Appendix A. 
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4-1- 2-loop results without resummation 

The 2-loop diagrams are shown in fig. [TR where we have ignored purely scalar diagrams 



since they are suppressed by A ~ e 3 and are lower-order than e 4 T 4 . Ignoring resummation 
for now, we shall discuss how to compute the diagrams. Fig. |15|a contributes 



^ V }») = -ieW yf 4Q 2 -4(P-Q) 2 /P 2 

M i 2 ^ %% (pa + M 2 )(Q 2 + m 2 )[( p + Q) 2 + m 2]_ l^.bj 

This may be rewritten in terms of integrals similar to I(m) and H(m) of the last section 
by repeatedly expanding factors in the numerator as sums of denominators, e.g. 

2(P • Q) = [(P + Q) 2 + m 2 ] - (P 2 + M 2 ) - (Q 2 + m\) + (M 2 + m 2 - ml), (4.7) 

and by noting that 



1 1 1 



P 2 (P 2 + M 2 ) M 2 P 2 M 2 (P 2 + M 2 )' 
In this way, (|4.6|) may be rewritten as 



(4.1 



= ~£ e ^ J ( M )[ J ( TO i) + J (^ 2 )] - I(roi)I(m 2 ) 



+ (M 2 - 2m 2 - 2ml)H(mi, m 2 , M) 

^2 _ (4.9) 
+ 1 ^ 2 2J [/(AO - HO)] [I{mi) - /(ma)] 

f 777, ^ 777^1^ — — 

+ m 2 [iJ(TOl1 m2 ' M) ~ ^ (mi ' m2 ' ° )] 



where /(m) is (|3.11|) as before and 



H(mi, m 2 ,m 3 ) = ^ 4e ]£ X 7^3— 2777=^ — "277775 , ^2 i 5r ( 4 - 10 ) 

TpTq (P 2 + m()(Q 2 + m 2 ) [(P + Q) 2 + m 2 ] 

H is simply a generalization of the H(m) of the last section to the case of unequal masses. 

We now need a high-temperature expansion of H. It is easy to relate H to H(m) by 

considering the difference 

1 



P(toi, m.2, 7713) — H(m) = /j 



4e 



p7qK^ 2 + m 2 ) (Q 2 + m 2 )[(P + Q) 2 + m 2 } 

1 



(P 2 + m 2 ) (Q 2 + m 2 )[(P + Q) 2 + m 2 ] 

(4.11) 
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for any choice of masses small compared to T. The difference is more UV convergent than 
either alone, and the contributions to the difference involving heavy (po ^ or qo ^ 0) 
modes are suppressed by the small masses. The dominant contribution to the difference 
comes from the three-dimensional piece (po = Qo = 0) and is order T 2 . The simplest way 
to evaluate it is to switch to configuration space,* and set e to zero since the difference 
( f4.f 1|) is UV convergent in three dimensions. The three-dimensional propagator is simply 
e~ mr j 'Anr , and we obtain 
H(mi, m2, m^) — H{m) 

T 2 f°° dr 



(4tt) 2 
T 2 



f°° dr 

/ — {exp[— (mi + m 2 + m^r] — exp(— 3mr)} + 0{mf 
Jo r 

{— ln(mi + m,2 + ms) + ln(3m)} + 0(m 2 ). 



(4.12) 



Adding this to the result ( |3.17|) for H{m) gives 

r Tl i m i +m 2 +m 3 ' 

H(mi, m 2 , m 3 ) = H 



1 T i(l + ii + l J?) +iln ( E 1 + 2 



647T 2 \e \T 2 J \mi + rri2 + 1713 

+ 0(m 2 )+0(eT 



(4.13) 



1 + Le +ln(^\ +121n^ ] + 6 + 4c B -3c H 



may now easily be expanded: 

a (a) = e 2 MT 3 _ e 2 M 2 T 2 
^ i 12(4tt) 24(4tt) 2 

+ 0(e 9 / 2 T 4 ). 

(4.14) 

The 0(e 9 / 2 T 4 ) comes from terms of order e 2 mMT. Though one could keep track of such 
terms in a 2-loop calculation, and so improve the error to 0(e 5 T 4 ), we have not bothered 
to do so. We shall focus only on the leading correction to the one-loop potential for 

e 4 < A < e 2 . 

The other graphs may be computed in similar fashion, and the results are given in 
Appendix A. 



* We thank Lowell Brown for this observation. 
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4-2. Resummation: effective masses 

We need to resum the propagators to include the dominant thermal corrections to 
the masses. We shall use the prescription of section |3.4j (method II) of only resumming 
the static fco = modes. In the Abelian Higgs model, the longitudinal (A$) mass at 
(po = 0, p — > 0) becomes 



+ 1st- 



(4.15) 



while the transverse mass remains the same, which we shall continue to denote by M: 



M 2 = eV- 



(4.16) 



The leading contribution to the scalar thermal mass changes the mf of 



rh 2 (<p) = m 2 (<j ) )+[h + h + 3eA^. 



to 



(4.17) 



The resummation of the 1-loop potential 



= - 
2 



gives 
1 



1 



3M 
64^ 



4 r 



In 



A 2 



2c, 



. —(2M 3 + Ml)T + 
+ 0(e 9 / 2 T 4 ) 



(4.18) 



T 2 ) 3 

We can now proceed as we did in section Before doing so, we should dwell a little 
more on the convergence of the loop expansion. (The following discussion is important 
for providing a well-behaved loop expansion but will turn out to have little practical 
importance for the actual computation of the 2-loop potential to (9(e 4 T 4 ).) Consider the 
multi-loop diagram of fig. ffBl The cost of adding loop A to the diagram is e • <p 2 T/M 3 ; the 



e (j) comes from the vertices, the T because it's a one-loop integral that's not quadratically 
divergent, and the 1/M 3 to make a dimensionless quantity. More specifically, the 1/M 3 
arises because adding the loop added three propagators to the diagram (two vector and 
one scalar) and a loop integral d 3 A; which is dominated by its infrared behavior. Both the 
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new loop A and the loop it is attached to are dominated by momenta of order M rather 
than m, and so each propagator gives 1/M 2 and the new d 3 k gives M 3 . The total loop 
cost of e 4 (j) 2 T/M 3 is the same order as e 2 T/M ~ e (see Table 1), which is the vector loop 
expansion parameter that we identified in the introduction. 

Ignoring the dominant e 2 T 2 piece that is absorbed by resummation, the cost of loop 
B is e 2 Tm/M 2 ~ e 3 / 2 since the loops are dominated by momenta of order m for the scalar 
and M for the vector. Similarly, loop C cost XT jm ~ e 3 / 2 after the dominant XT 2 piece 
is absorbed by resummation. Adding line D cost e 2 T/M ~ e. Adding line E, however, 
cost e 2 TM/m 2 ~ 1. So it would seem that the loop expansion parameter is not order 
e but instead order 1. The problem arises because of the difference of scales between m 
and M. When we defined the effective scalar masses rh in ( [4.17| ), we only included the 
dominant 0(e 2 T 2 ) or 0(XT 2 ) contributions of vector or scalar loops to the self-energy. 
However, the sub-leading 0(e 2 MT) vector-loop contribution is the same size as rh 2 itself 
since rh 2 ~ e 3 T 2 ; it is not a perturbation and so must also be included in the resummation. 
(Yet higher-order terms are perturbations and so do not need to be included.) The new 
procedure for resummation may be viewed as follows: first integrate out all of the heavy 
modes to obtain an effective theory for momenta ko = and I; <T, next integrate out 
the vectors (and scalars with k ~ M) to obtain an effective theory for momenta ko = 
and k <C M, and only then finally integrate out scalar loops controlled by k ~ rh M. 
This procedure is equivalent to modifying our resummation prescription by including the 
0(e 2 MT) corrections to the one-loop scalar mass (^.17| ): 



This result may be computed directly from one-loop diagrams or, more simply, from the 




(4.19) 



second derivative through 0(e 3 T 2 ) of the resummed one-loop potential ( |4.18| ). Using the 



resummation ( |4.19|) , the cost of adding a new loop will now always be <e. 
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Diagrammatically, this resummation corresponds to the dominant pieces of a subset 
of "once-iterated" ring diagrams, as shown in fig. [H] for the improved scalar contribution 
to the one-loop potential.* The smallest loops of the diagram are hard, with momenta 
of order T, the next smallest are vector loops with momenta of order M ~ eT, and the 
large Higgs loop is softer yet, with momenta of order m ~ e 3 / 2 T. Because of the hierarchy 



of scales, it is a good approximation at each level of fig. [17] to approximate resummed 
propagators l/[p 2 + U.(p)] by l/[p 2 + 11(0)]. 

Though the foregoing discussion was necessary to establish an adequate procedure for 
obtaining a controlled loop expansion, the details of resumming the scalar masses turn out 
not to be relevant at the order under consideration. Though the results of some graphs 
have potentially significant terms of the form e 2 mT 3 , which would be affected by the details 
of how m is replaced by in, all such terms cancel in the final potential. The m 2 T 2 term 
of the one-loop potential is not modified to rh 2 T 2 if we only resum ko = modes (method 



II of section |3.4p.t All other terms involving m are lower order than 0(e 4 T 4 ). [We would 



need to take care to use the correct resummation ( [4.19|) for m if we were keeping track of 
0(e 9 / 2 T 4 ) contributions to the potential.] 



4-3. Resummation: 2-loop diagrams 

Now let us turn to the resummation of actual diagrams. The resummation of fig. |T5|a 



is depicted in fig. [18 



In the last pair of terms of fig. |18|b, resumming the vector line makes no difference 
because the Aq polarization couples to the ko components of the scalar momenta, which 
are zero. The only difference between the last two diagrams is therefore the scalar masses. 
It is easy to see by examining the general result (|4.9|) for the setting-sun diagram, and 



* The discussion of such diagrams appears in a related context in Ref. [26|. 
' If all modes are resummed (method I), the m 2 T 2 term is modified only by the addition of 
a constant, which we would ignore. 
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remembering that the three-dimensional version ( |3.36 ) of I(m) is simply Is(m) = — mT/47r, 
that these three-dimensional diagrams do not depend on the scalar mass at 0(e 4 T 4 ). So 
their difference is ignorable and may be dropped. 

In the other two pairs of terms of fig. |18|b, the heavy loops simply act as 0(e 2 T 2 ) 
mass insertions in the light loops. For example, the first pair of terms contributes to the 
potential 



--e 2 T 
2 



-e 2 T 
2 



d 3 p 



2tt) 3 [p 2 + M 2 p 2 + M 2 _ 

d 3 p 
(27F 

d 3 p 



p 2 + M 2 p 2 + M 2 



90^0 

qo^o 



d 3 q (P + 2Q) lI (P + 2Q)„ 
[2tt) 3 {Q 2 + m 2 )[(P + Q) 2 + m 2 ] 

d 3 q 4Q^Q„ 
(2n) 3 Q 4 



[l + 0(p 2 /Q 2 ,m 2 /Q 2 ) + ---] 



(2nY 



p 2 + M 2 p 2 + M 2 



x _i e 2 T 2 [i + (p 2 /T 2 ,m 2 /T 2 ) + ...} 



1 



e 2 (M L - M)T 3 + 0(e 5 T 4 ' 



12- 



(4.20) 



where n M = (1,0) and P = (0,p). It is the MT 3 and mT 3 terms which resummation is 
supposed to eliminate. This does not occur on a graph-by-graph basis, but they cancel in 
the sum over graphs. The M^T 3 and mT 3 terms cancel against the counter-term graphs 



of fig. 15ti. Evaluating the rest of fig. lSb similarly, one finds 



1 



2ey (a) = 2ey (a) _ 
P ^ $ 12-47T 



e 2 (M L - M)T 3 + 0(e 5 T 4 ) 



(4.21) 



The results for the other diagrams may all be found in Appendix A. When added 
together, they give the full, improved 2-loop potential: 



^ 2 (T)Lu 2 (T)+(h(T) + 3e 2 (T)\ ^ 



4rp2 r 



e 4 T : 



\2M h J \2M 



+ 
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±-{2M 3 + M 3 )T+^{T) 



A(T) + 



9 

36e 4 
(4tt) 2 



(4.22) 



c B + - 



+ 0(e 9 / 2 T 4 ) 
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where 



e 2 (T) = e 2 + 




(4.23) 



A(T) = A + 




(4.24) 



0(T) = 1 + 



is 2 + 0(e 5 T 2 ). 




(4.25) 



(4.26) 



It is important to remember that this result was derived in the high-temperature limit, 
and the T — > limit of the terms shown in Q4.22| ) is not the same as the zero-temperature 
potential. The T — > limit of ( |4.22|) is trivial since the \(T)<p 4 term vanishes if the scale 
of A is run to zero. For the true potential at zero temperature, the relevant coupling is 
instead \(<p). This distinction is important to keep in mind when using the 2-loop result 
for numerical work. 

4-4- The vacuum shift 

How much do the 2-loop corrections just computed affect the asymmetric vacuum at 
the phase transition? The simplest quantity to examine is the vacuum expectation value. 
This VEV is not a physical, gauge-invariant quantity, but the size of its shift due to 
2-loop corrections will still give us a good test of whether the loop corrections are small. 
In a companion workjH one of us has instead computed the magnetic screening mass in 
the asymmetric vacuum, which is physical and gauge invariant. (Actually, we suspect that 
the VEV may be gauge-invariant to the order we have computed, but we do not know for 
sure.) 

Normalize the potential so that V(4>, T) is zero at the origin. The VEV </> of the 
asymmetric vacuum at the critical temperature T is then determined by simultaneously 
solving 



V(cf>,T) = 0, 



d 4> V{<(>,T)=0. 



(4.27) 
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Write V = V\ + AV where V is the full potential, V\ is some approximation (say, the 
one-loop potential), and Ay is a small perturbation. Write = 0i + A0 and T = T\ + AT 
where 4>± is the asymmetric vacuum of the approximation V\ at its critical temperature 
Ti; that is, 

Vi(0i,Ti) = O, fyVi(0i,Ti) = 0. (4.28) 
Linearizing the equations (|4.27|) in A(p and AT, one finds 

^ AVd.dTVx - (d,AV)(d TVl ) 
^ (d^idrVx) ' (4 -' Uj 

where all the quantities are evaluated at the VEV. Note that A</>/0 is order AV/V, where 

V is the typical size of terms in the potential, and so A<p/(f) is of order the loop expansion 

parameter. 

Now identify V\ with the one-loop potential and note that 

d T V x = \e 2 T$ 2 + 0(e 3 T 2 ). (4.31) 
Plugging into the solution Q4.30|) for A<p: 



¥ = -4^(^)+0(e^). (4.32) 



It is important when applying this formula to remember that AV((p : T) should be replaced 
by AV(4>, T) — Al/ (0, T) if the potential is not already normalized to zero at <p = 0. This 
result has the important property that is vanishes at leading order if Ay is proportional 
to (j) 2 . So the only two-loop contributions to the potential ( |4 . 2 2| ) which contribute to the 
(leading-order) shift in the VEV are those that involve logs of masses: the e 4 T 2 2 lnM(0) 
and e T 2 (p 2 InM-^^) terms. (The e 4 4 terms, which would also contribute if taken as part 
of Ay, came from the one-loop contributions.) Equation (|4.32[) applies to all models we 
shall examine. 
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For the Abelian Higgs model, ( |4.32| ) reduces to 



e 4 T 2 / M 2 \ 

2 + 772 > ( 4 - 33 ) 



(4yr) 2 m 2 V M 2 

where Ml, M and fh\ are as usual the effective masses at the phase transition given by 
( [11 5|) , (|4.16| ) and (f4.19| ) , and where the one- loop approximation to the VEV <fi is sufficient 
for the right-hand side above. The one-loop approximation to <fi doesn't have a simple form 
in terms of T, M W (T = 0), and m^T = 0) [at least not for the whole range e 4 < A <C e 2 ], 
and is best computed numerically. 



5. Fermions 



Fermions are simpler to deal with than bosons because they do not require resum- 
mation. This is because Euclidean fermions have frequencies ko = (2n + 1)itT which are 
never zero, and so their self-energies may always be treated perturbatively. For the same 
reason, the two-loop contributions to the potential that involve fermions will never have a 
logarithmic dependence on the fermion mass mf (</>) and so will not affect the shift in the 
VEV. Two-loop fermionic diagrams do contribute to the effective potential at 0(g 4 T ), 
however, and so we shall treat them here (and later in our discussion of the Minimal 
Standard Model) for completeness. 

As a simple example of a theory analogous to the weak interactions, let us chirally 
couple a single Dirac fermion to the Abelian Higgs model: 



— ie/i 

— iefi 



1-75 
2 

1-75 



l-7£ 



1 + 75 



(5.1) 



We shall treat the Yukawa coupling g y as being of the same order as the gauge coupling 
e. The fermion mass in the presence of a background field <p is 

1 



rrif = 



V2 
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(5.2) 



and the one-loop potential now picks up the familiar fermionic contribution: 



fi 2e V (1) -> ^ 2e F (1) + 4J f (m f ), (5.3) 

where O 



i _^ f > 

i 



Ji{rn i ) = ~^ 2 % ln(K 2 + m f 2 ) 



1 o„o I /I fl , , / \l 2 



€ 



- + In ^ - 2c F 



2^2 



(5.4) 



= const. + -m f 2 T 2 + ^ 

+ 0(m 6 /^ 2 )+0(e), 
and where we have defined the constant 

c F = In 7T — 7 E = c B — 2 In 2. (5.5) 

The superscript (f) on the integral indicates that the frequency sums are over fermionic 
momenta ko = (2n + 1)tcT. The factor of four in ( |5.3|) corresponds to the four components 
of a Dirac fermion. As is conventional in MS regularization, we take the Dirac trace of one 
to be tr 1 = 4 rather than 2 2_e . 

Our working definition of 75 in dimensional regularization is the naive one that@ 

{75,7m) = 0, (5.6) 

7l = I- (5-7) 

This definition is adequate for calculations that do not require traces involving an odd 
number of 75's. Additional properties are 

tr(757M7^7o-7r) is anti-symmetric; (5.8) 
tr(757M7^7<r7r) = 4ie MI/0 . T + 0(e) x ambiguity 

on the four-dimensional subspace fi, v, a, r = 0,1,2,3. (5.9) 

The anomaly is intimately related to the fact that the ambiguous term above cannot be 
explicitly defined. In our calculation, however, this term will never be relevant. 
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The new 2-loop diagrams are shown in fig. [H| As an example, take fig. 19 i and expand 
numerators in terms of denominators as we did in ( |4.6| ) through (^^) for the Abelian Higgs 
model: 

2nr (U = 1 r 2^f (f) ^f (f) [S»u ~(P + QUP + QU(P + Q) 2 ] 

^ i 2 Tp Tq (P 2 + m f 2 )(Q 2 + m 2 )[(P + Q) 2 + M2] 
x tr(-iP + mt ) 7/ x (^-^J (W + mfh* P ~ 

(5.10) 



(f) ' (f) [^-(P + Q) M (P + Q),/(P + Q) 2 ] 



8 5q (P 2 + m f 2 )(Q 2 + m 2 )[(P + Q) 2 + M 2 ] 

= e 2 {(l - e)[/ f (m f )] 2 - 2(1 - e)/ f (m f )/(M) 

+ [(1 - 2e)m 2 - (1 - e)M 2 ]# ff (m f , m f ; M)}. 



If and are defined analogously to the bosonic I and H of previous sections: 



/•(f) 



1 



(5.11) 



P/ ff (m f2 ,mfi;m) = // 



4e 



h K %{K* + m f 2 2 )(Q 2 + m 2 1 )[(P + Q) 2 + m 2 ] ' 
The high-temperature expansion of Jf(raf) may be obtained from ( |5.4| ) using If 



(5.12) 



-m f 1 d J/dmf : 



/f ( m )=_l T 2 - -^ 
u ; 24 16tt 2 



m 



2c, 



+ 0(m 4 /T 2 ) 



(5.13) 



In our earlier calculations, it turned out that the coefficient L e of the eT 2 term of I(m) 
was unimportant because it canceled in the final answer. When fermions are added to the 
theory, one finds that the final result does depend on the combination — t e , which it 



behooves us to calculate. One may find by a derivation of If similar to that for / in ( |3.32j ) 
that 

4 f) =t £ -21n2. (5.14) 
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ff) 

We shall always make this substitution for L e , and in this way the 6 e 's will cancel in the 
final answer for the effective potential. 

The leading 0{T 2 ) contribution to turns out to vanish, and 



Hff(mf 2 , run; ™) = 0(mT) + 0(eT 2 ). 



(5.15) 



A derivation is given in Appendix B. The consequence of (|5.15|) is that the term in 
( |5.10| ) may be ignored when studying the effective potential through 0(e 4 T 4 ). terms 
are similarly ignorable in all diagrams we shall calculate. 

The high-temperature expansion of our result ( |5.10| ) for fig. |19]i is then 



e 2 MT 3 
12-4% 
e 2 M 2 T 2 



+ 



e 2 m 2 T 2 



4(4tt) 



7 + te + ln (^ 



7 + ^ + ln (^2 



1 - 2c P 



l-2c P 



21n2 



-ln2 
3 

+ 0(e 5 T 4 ) 



(5.16) 



12(4tt) 2 

The results for the rest of the diagrams, and the results of resumming the ko = modes 
of bosonic propagators, are given in Appendix A. The final 2-loop potential is then 
V {2) = [Abelian Higgs result of ([122])] 



+ 



l 

12 

e 2 9 2 y T 2 



(47T) 



1 1 , 

-c F ln2 

2 F 3 



+ 4! ( 



(4tt) 2 
(4tt) 2 



1 

12 
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1 , \ e 4 T 2 
2 CF -6 ln2 J + (47F 



11 1, 

h -c B + - In 2 

18 3 3 



(5.17) 



where the running couplings A(T) and so forth in (|4.22| ) now include fermionic effects. (See 
Appendix A for formulas for the running couplings.) 



6. Nonabelian Theories 

Extending our analysis to nonablelian theories is fairly straightforward; it merely 
involves the computation of some new graphs, shown in fig. |20|. In Appendix A, we give 
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results for the SU(2) Higgs model with a single Higgs doublet, defined by 



C = --F^ u F a " u + |L>$| 2 - 
where $ is an SU(2) doublet and 



F(|$| 2 ), V(\$\ 



11 3! 1 



du-i^gAp-r } <l>. 



(6.1) 



(6.2) 



We once again normalize (p to (p 2 = 2|$| 2 . 

Reducing the graphs of fig. |2(] to simple scalar integrals like I(m) and H(mi, m 2 , ms 
turns out to require the introduction of a new function: 

{P-Q? 



L(mi, m 2 ) 



- p ^P 2 (P 2 + m 2 )Q 2 (Q 2 + 



mi 



(6.3) 



One can derive a high-temperature expansion of this function, which we give in Appendix 
B, but it is unnecessary because terms involving L(mi, 1712) cancel exactly between figs. |2( 



and |20jo. 

The only part of the calculation of the two-loop potential that is not a completely 
straightforward extension of our previous model calculations is the resummation of fig. pOm. 



The resummation is shown in fig. |2l|. Recall that only the longitudinal (Aq) polarization is 
resummed, and split propagators into longitudinal (L) and transverse (T) parts. By using 
the fact that the three-vector coupling vanishes for LLL and LTT when all frequencies po 



are zero, fig. ^lj may be rewritten as fig. ^2|. The three dimensional diagram in the last 
pair of terms is easy to evaluate, and the result is given in Appendix A. 
The total ring-improved 2-loop potential is given by 
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where we have introduced the notation 



M* = Ml - M 2 (6.5) 



for the ^-independent thermal mass. Note that though the M-^M term above is linear in 
4> as <fi — > 0, this linear term cancels against the <fi — > behavior of the M 2 cbyc ln(2ML +M) 
term. 

We can now compute the relative shift in the VEV from its one- loop value using Q4.32Q : 



A0_, g*T 2 [87 , 3M 2 byc 3M 2 3M£ bye , ( 2M h + M 



(j) (An) 2 m{ 



' clcbyc acbyc 



(6.6) 



32 2M L M 16M 2 M 2 V 2M debye 

where M, Ml, and mi may be evaluated at the phase transition using the one-loop ap- 
proximation. rh\ is obtained by taking the second derivative of the one-loop potential 
through 0(g 3 T 2 ): 

~ 2 2 (, 9 2 \ T 2 3g 2 MT 3g 2 T ( M 2 \ , w , 

™? = mf + (a + i9 2 j T2 - - v- - w U: + Ml ) • (6 ' 7) 

Note that the MlM term in the potential contributes to A<p/<p, whereas only logarithms 
contributed in the Abelian case. 



7. The Minimal Standard Model 

Finally, consider the Minimal Standard Model with a single Higgs doublet, which 
differs from the previous models only by the complexity of keeping track of the various 
couplings. Our conventions for coupling constants are such that 

D» = a M + ±g 2 A„ ■ r + l -Y gi B^ (7.1) 

the hyper charge is normalized so that Q = T3 + Y/2, and the Yukawa coupling is 

9 y Ql ■ $ *r + h.c. (7.2) 
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where $ is the full complex doublet. g y is the top quark coupling, which is the only 
Yukawa coupling that we treat as non-zero. We shall also include QCD interactions for 
the quarks, with coupling g s . The couplings gi, g%, g y , and g s are formally treated as all 
having the same order of magnitude g. nf represents the number of families and is three 
in the Minimal Standard Model. 

At zero temperature, the mass matrix for the Z boson and the photon, in the (A^ 3 \ B) 
basis, is 

M 2 = ( \ 9 ^ 2 -J9i92<p 2 \ _ f cos9 w sin 9^\ f M 2 \ /cos# w -sin# w \ 
V — i^i^^ 2 \gl4> 2 ) \-sin0 w cos9 w J \ J \sin9 w cosd w ) ' 

(7.3) 

When we include the Debye masses generated for the longitudinal (Aq) components at 
finite temperature however, one finds© 



L "V -bi92^ 2 \gl<P 2 + {\ + *-t)glT 2 

cos 9 sin 9 \ ( M^ z \ ( cos 9 — sin 9 



- sin 9 cos 9 ) \ ) \ sin 9 cos 9 



(7.4) 



and we have diagonalized to define an effective mixing angle 9(<p,T) and effective longitu- 
dinal masses M^ z ((f),T) and Ml 7 (<^>, T). 

The result for the two-loop potential is given in Appendix A, where we list the contri- 
bution from each of the graphs of fig. [23| Only those diagrams which give a ^-dependent 
contribution to the potential at 0(g T 4 ) are shown. Note that there is a QCD correction 



in fig. ^3|i which gives a potentially large contribution to the potential of order g 2 g 2 (p 2 T 2 



but which, like all fermionic contributions, does not contribute to the leading-order shift 
in the VEV from its one-loop value. 

It is cumbersome to write down the total answer for the two-loop potential, and we 
shall leave it in the form of Appendix A where each diagram is given separately. The 
reader should note that the relative shift in the VEV is given by ( |6.6| ) in the limit that 
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<7i — > 0. This is because, as discussed earlier, the fermion contributions do not affect the 
shift in the VEV except indirectly through their effect on the one- loop values of M, Ml, 
and rhi at the phase transition. For example, rh\ in this limit changes from (|6.7| ) to 



2 2 (, 9 2 n ,\ T 2 3g 2 MT 3g 2 T ( M 2 w , 



8. Numerical Results 



As a test of the size of loop corrections, we shall now compute the shift A<fi in the 
VEV at the phase transition due to the inclusion of two-loop corrections in the Minimal 
Standard Model. We do this numerically, by running couplings to Ji = T and then finding 
the phase transition for both the one-loop and two-loop approximations to the potential 
given in Appendix A. 



As discussed in section |4.4| , the relative size of this shift is O(g), the loop 

expansion parameter. This means that we need to specify the original zero-temperature 
parameters of the theory, such as A and g 2 , to relative order g in terms of physical quantities 
such as the zero-temperature Higgs and vector masses. Tree-level relations are adequate 
for most quantities, since corrections at zero-temperature are suppressed by g 2 instead of 
g. The exceptions are A and v 2 since they are 0(g 3 ) by our power-counting convention 
and receive 0(g 4 ) corrections at one-loop. To obtain A and v 2 in terms of physical masses, 
consider the one-loop potential at zero-temperature: 
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We shall refer to the zero-temperature VEV as a. The physical Higgs mass at the order 
under consideration is just the second derivative of this potential at a* Solving for A and 
v 2 in terms of the physical masses, one finds 
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3ml 



3 [ 3 4 
32tt 2 V 2 
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2 h 64tt 2 
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1.2) 



1.3) 



where the bars denote that masses are evaluated at the VEV a at zero temperature with 
renormalization scale p,, and where a may be expressed as a = M w /g. 

Our numerical results for the shift A</>/</> due to 2-loop corrections are shown in fig. 24 
as a function of mh, assuming a top quark mass of 100 GeV. To get the solid line, we 
computed cf> and T c independently for both the one-loop potential and for our full result 
for the 2-loop potential. To help control errors of the one-loop computation at very small 
mh, we have used exact results for the one- loop fermion and vector contributions rather 
than high-temperature expansions. 

Alternatively, the dashed line gives the simpler calculation of the pure SU(2) result 
( |6.6|) for the shift, which we have evaluated using the Minimal Standard Model results for 
M w and Ml w in place of M and Ml. We have also replaced fh\ by the complete second- 
derivative of the one-loop potential rather than by ( |7.5| ) because, due to fine cancellations 
among terms, the assumption that g\ — > used in (|7.5| ) turns out to be a very bad 
approximation if the one-loop VEV and critical temperature were computed with gi ^ 0. 



* The physical Higgs mass is actually determined by the pole of the Higgs propagator and so 
is given by the solution to the real-time dispersion relation P 2 — M 2 = U(P 2 ) where n is the self- 
energy. The claim that it is given by the second derivative of the effective potential corresponds 
to the approximation H(P 2 ) — > n(0) in this dispersion relation. The difference H(P 2 ) — n(0) is 
order g 5 a 2 and does not affect our derivation of A through 0(g 4 ). 
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The dashed line should approximately match the solid line when both (1) is small 

and (2) is not so small that the high-temperature expansion has broken down. This 
correspondence is evident in fig. E3. The deviation of the lines at large is a manifestation 



of the breakdown of the loop expansion. 

The dotted line in fig. |2| is a diagnostic of our ubiquitous assumption (A ^> g 4 ) that 



we may use the high-temperature expansion. It is the ratio of the 0(M 6 /T 2 ) piece of the 
one-loop potential to the 0(M 4 ) piece. This ratio is taken in the asymmetric vacuum at 
the phase transition. More specifically, 

v * = ttt^ (6c B M 4 + 3c B M z 4 - 12c F m 4 t ) , 

C(S) (8 ' 4) 

V6 = 768^ (6M ' + 3M ' " 7 ' 12Wt6) ' 



and the dotted line is l^e/^l- One can see in fig. £4] that the high-temperature expansion 
begins to break down at small m^. 



Fig. ^5] shows the result for <j)/T at the phase transition, which is the quantity used in 



Ref. |E| to extract the upper bound on the Higgs mass for weak baryogenesis. According to 



Ref. [12 1, a necessary requirement for a successful scenario of baryogenesis is that <f>/T > 1. 
The dashed and solid lines show our one-loop and two-loop results respectively. At the 
experimental lower bound of 60 GeV, the one-loop result of 0.47 is inadequate. The 2-loop 
corrections boost <j>/T by 40% to roughly 0.65. This remains inadequate. Unfortunately, 
the corrections are large enough that this conclusion does not impress us as airtight — the 
validity of the loop expansion is only marginal here. Note that, for a consistent comparison 
at this order, the sphaleron mass and resulting limit on <p/T should also be consistently 
determined to the same order. 



Figs. ^6| and fig. ^7] show our results for a top mass of 180 GeV. In parting, we 
remind the reader that the Minimal Standard Model is only a specific testing ground for 
these issues and that the particular constraints just discussed are evaded in multiple- Higgs 
models. 
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Appendix A. Results for 2-Loop diagrams 

In this appendix we collect 2-loop results for individual diagrams in the various theories 
discussed in this paper. Since the theories involve many of the same diagrams, and differ 
only by group factors, we only compute them once. The D functions defined at the end 
of this appendix are the results of diagrams where we have factored out a combination 
of coupling constants, symmetry factors, and group factors. The subscripts of the P's 
indicate the types of particles involved in the diagram: S for scalar, V for vector, F for 
fermion, and r\ for ghost. Note that that the effective potential and the sum over diagrams 
have a relative minus sign in Euclidean space. We shall display only the (^-dependent 
pieces of contributions to the potential and shall not explicitly indicate the presence of 
temperature-dependent constants by "const." For brevity, we shall drop the unimportant 
/x 2e terms that keep track of the dimensions of the potential V and the field <p; but every 
instance of V and <fr in this appendix should be understood as /j> 2e V and Throughout 
this appendix we adopt the notation 



1 



(1 + eie), 



(A.l) 



e 



where t e is the unimportant constant of ( |3.12| ), and we define 



Hs(m 1 ,m 2 ,m 3 ) 



1 



T 2 ln(mi +m 2 + m 3 ), 



(A.2) 



(4tt) 2 



to be the (non-constant) three-dimensional piece of H, as in ( |4.12| ). 
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To actually use the formulas of this appendix for computing two-loop corrections 
numerically, sum the individual contributions but (1) drop all 1/e terms, since they cancel; 
(2) set ft to T and so drop all ln(/i 2 /T 2 ) terms; and so (3) replace couplings g, masses v 
and fields <p everywhere by g(T), v(T) and 4>(T). 

Unresummed results for the diagrams in (A.l) and (A. 2) may be obtained by replacing 
all £>'s by V^s and dropping V^ h K 

A.l. Abelian Higgs Model 
Counter terms: 

</>bare = (l + 2 (4^ e H ) ^ bare = ( 1 ~~ 6(47r)2 e ^ ) A 

2 _ ( 1 i e 2 , "\ „2,,2e \ _ ( \ , 18e 4 , "\ ,,2e 



^are = (l + + ' " j e A bare = (A + ^ + • • • ) ^ (A.3) 

= (l " + " " ) ^ 

Multiplicative corrections are shown only through order e 2 , and order A corrections are left 
out since A ~ e 3 by assumption. Running couplings and fields may be found by replacing 
1/e -> ln(T 2 /fi 2 ) and </> bare -> 0(T), e 2 are -> e 2 (T), etc. above. 

Masses: 

M 2 = e 2 ^ 2 , m\ = -v 2 + -\<f> 2 , m 2 = -v 2 + -A0 2 (A.4) 

2 (3 

M 2 = e 2 <p 2 + le 2 T 2 , rh 2 = m 2 + Qa + 3e 2 ) ^ + ■ • • (A.5) 

We shall not bother showing the sub-leading corrections to the thermal scalar mass, as in 
( f4.19|) , since they cancel out in the final result for the potential to 0(e 4 T ). 

1-loop result: See eq. ( |4.18| ). 



2-loop diagrams of fig. 15: 



V (a) = _ I e 2 Pssv(TOl7 ^ M)) ( A . 6a ) 

V {b) = -^e 4 (/) 2 ^svv(mi, M, M) (A.66) 
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= -^e 2 [I?sv(TOi, M) + V sy (m 2 , M)] 


(A.6c) 




_ e 2 M 2 T 2 n_2\ 5y4 
24(4tt) 2 Ve 3; ^ J 


(A.6d) 




8(4tt) 2 Ve 3 / 


(A.6e) 


yd) 


= 0(e 5 T 4 ) 


(A.6/) 




= eVT 2 l 11/2 4 
2(4tt) 2 e + 1 J 


(A%) 




e 2 M L T 3 , e 2 (m 1 + m 2 )T 3 , ^ 9/2 ^ 
" 24tt ' 32tt + 1 J 


(A.6/1) 



A. 2. Chiral Abelian Higgs Model 



Counter terms: 



'bare 



bare 



"bare 



bare 



1 + 



(3e 2 -^) 



2(47r) 2 e 



+ 



4(4tt) 2 



1 + 



+ • 



(47r) a e 
(3e 2 - ff ;) 

(4tt) 2 £ 



+ 



e 2 ^ 



^4-bare 




V'R.bare 




Abare 




fi'yjbare 


= ri + 



2(4tt) 

2 
fv 



+ ■ 



4(4-7r) 2 e 
(18e 4 -6g4) 



(47r) 2 e 



(29 2 -3e 2 ) 



2(47r) 2 e 



+ 



2e 



(A.7) 



Formulas for running couplings may be obtained as in section ~KA 



Masses are as ( |A.4| ) with: 



nit 



2 ±2 



Mi 



,2 ±2 



— i— - e 2 T 2 
2 



m 2 



m, 2 + 



jjA + 3e 2 + s 2 )^ + ---. (A.8) 



1-loop result: 



yd) = 



( -A + 3e z +^ 



3M 4 



4(4tt) 2 
+ 0(e 9 / 2 T 



In 

4\ 



J"2 



2 
3 



1-1 


1 


12 


^ 127T 




"1 4 

+ (4tt) 2 


-2c B 



(2M 3 + M[)T + ^X(f> 4 



(A.9) 



2-loop diagrams of figs. |i~5|a-c same as (A.6.) 
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2-loop diagrams of figs. 



pld-h: 



e 2 M 2 T 2 (\ 2 



+ 0(e 5 T 4 " 



8(4tt) 2 Ve 3 
e 2 (3e 2 -g 2 )cP 2 T 2 fl 2 



8(4tt) 2 
= 0(e 5 T 4 ) 

(3e 4 - g^ 2 T 2 1 



§-§)+0(e 5 T 4 ) 



6(4tt) 



+ 0(e 11 / 2 T 4 



e^ + (3e 2 ^ 9 2 ) {rni+rn2)T3 + o{e9/2T4] 



16n 



96tt 



(A.10d) 
(A.fOe) 

(A.10/) 

(A.lOsr) 
(A.10/i) 



2-loop diagrams of fig. f9: 



V® 



^e 2 "D F Fv(wf,mf,M). 



1 



1 



7# y £ > FFs(«Tf, mf,mi) - ^fi'yl'FFS^f, m f ,m 2 ) 



12(4tt) 



Q -2 In 2^ + 0(e 6 T 4 ) 



V« = 



(A.IOz) 

(A.lOi) 
(A.lOfc) 

(A.10Z) 



A3. Gauged SU(2) Higgs Theory 



Take # s = g y = g x = e = 6> w = 9 = 0; set M z = M w = M, M Lz = M Lw = M L , # 2 = and 

ignore Ml 7 in the Minimal Standard Model results below. 
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A. 4- Minimal Standard Model 



Counter terms: 



Y'bare 

^4-bare 
-^bare 

^L,bare 



92, bare 
2 

^l,bare 
fi'yjbare 
Abare 



V, 



bare 



1 + 
1 + 



(9g| + 3gj - 12g$) 
8(47r) 2 e 

25 2n{\ g\ 



12 3 J (47r) 2 e 
1 10n f \ g\ 



+ ■ ■ 



+ 



A 



' ' 12 + 9 ) (47r) 2 e + ' 



B 



1 



4, 



4(47r) 2 e 

,2 



+ • • • Vl 



^R,bare — 1 



+ 



2(4vr) 2 e 
43 4n f 



R 



43 



6 3 7 (47r) 2 e 
1 20n f \ <? 2 



1+, 6 + 9 



) (4 



7r) 2 e 



+ 



+ 



2 2e 
02^ 



2 2e 

9iH 



9 
4' 
27 



1 + ( l9y ~ ^9s - «02 - ~—*9t 



17 
24' 



(47r) 2 e 



+ 



I -I (^2 4 + ^ + ^i-18^ 4 



(9^ 2 + 3g 2 - Ylg] 
4(47r) 2 e 



H v 



y ' (47r) 2 e 
2 



+ 



A// 



2e 



(A.11) 



Formulas for running couplings may be obtained as in section [A.l . 



Masses: 



M 2 = \gW, 


M z 2 = ^2 2 + S?)</> 2 , m f 2 = ^> 2 


(A.12) 


2 2 
r?^ = —v 


+ ^A0 2 , = -v 2 + ^A</> 2 


(A.13) 


ML 




(A.14) 


~ 2 2 , 
m i = m, + 




(A.15) 
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\gW + (| + f ) g\T^ 



W + (J + ^ 0?T 5 



4 

cos6>_ sin6> 
- sin 9 cos 9 



ML 



Ml, 



cos 9 
sin# 



— sin# 
cos 9 



1-loop result 

yd) = _ 



3 

V 



T21 

12 



12tt 



{AMI + 2M Z 3 + 2M L 3 w + ML + M^)T + -X^ 



3(2M 4 + M z 4 ) 
4(4tt) 2 

+ OG7 9 / 2 T 4 ) 



+ 



3m 4 
(47)2 



2c, 



(A.16) 



(A.17) 



2-loop diagrams of fig. |23|: 

cos6> w is defined below by the tree-level relation tan6> w = gi/g%, and e is <?2sin6> w . We 
shall also use the short-hand notation 

AM 2 = 2M Lw + M Lz cos 2 9 + M Ll sin 2 9 - 2M W - M z cos 2 W , 

(A.18) 

AM i = M Lz sin 2 + M hl cos 2 - M z sin 2 W . 

Take care to note that many of the graphs below are written in terms of the instead 
of the V. This is because of notational complications caused by the difference between 6> w 
and 9. 



y {a) = -i^ssvK» TO 2,M w ) - ^ 2 D4v(™2,m 2 ,M v 



8 cos 2 ft 



8 cos 2 e x 



■(1 -2sm A 9 w yV f ssv {m2,m 2 ,M z ) 



+ 



~e 2 2^ sv (m 2 , ma, 0) - ^-(g 2 2 AM 2 + g 2 1 AM 1 ) + 0(g 9 / 2 T 4 ) 
' 2 ^^svv(mi,M w ,M w ) 

I 4 / 

64 ^h~ [V * vvimi ' Mz ' Mz) " 4 ^ 3(mi ' Mz ' Mz) ] 



+ ^(7 2 e 2 tan 2 w [££ vv (m 2 , M w , M z ) - 4tf 3 (m 2 , M w , M z 



(A.19a) 
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1 

1 



+" ^le 2 [^vv( w 2,M w ,0) -4i? 3 (m 2 ,M w ,0)] 



+ —(g 2 cos 6 + gi sin0) 4 # 3 (mi, M Lz , M Lz ) 
+ - (g 2 cos + #i sin 0) 2 (g 2 sin - 01 cos 0) 2 H 3 (mi , M Lz , M Lt ) 

o 

+ ^(02 sin0 - 0i cos 0) 4 i? 3 (mi, M Lt , M Lt ) 

+ ^0 2 2 2 sin 2 0i? 3 (m 2 , M Lw , M Lz ) + ^g\g\ cos 2 0tf 3 (m 2 , M Lw , M Lt ) 
+ O(0 5 T 4 ) 

^ (C) = -^ 2 2 ^ V (^2,M W ) - ^[V^ v ( mi ,M w ) + vf v (m 2 ,M w )] 



(A.196) 



" Tfl-C 1 " 2 sin2 ^w) 2 ^ v (^2, M z ) - Je 2 ^ v (m 2 , 0) 

o COS t'w ^ 

~ h oSe^ [v ^ {mu Mz) + ^ v(m2 ' Mz)] ~ g^ 2 ^ 2 + 9 " AMl) 

- (192 + \g^j (mi - mi + 3m 2 - 3m 2 )^- + O(0 9 / 2 T 4 ) (A.19c) 



T 2 



4(4tt) 2 
+ O(0 5 T 4 ) 



(A.19d) 



F (f ) = O(0 5 T 4 ) 



9 4 i 3^ 2 ^ 2 i 3 t | 

r 



#2 + ^f2f 1 + gf 1 - 



T" 1 



8(4tt) 2 e 



+ O(0 9 / 2 T 4 ) 



(A.19/) 
(A.190) 



^ (h) = <7 2 2 Q + y) (2M Lw + cos 2 0M Lz + sin 2 0M Ll )^- 



1 5nf 



T 3 



+ 0f ( g + -g 1 ) (sin" 0M Lz + cos" 0M Lj ) — 



8tt 



+ (^f + ^i 2 + ^ 2 ) (mi + 3m 2 )|^ + O(0 9 / 2 T 4 ) 
y« = -40 2 ^ FV (m f ,m f ,O) 

[2^ FV (m f , m f , M z ) - p£ fv (0, 0, M z )] 



(A.19/i) 



3 2 2 
16 cos 2 V 



1 - - sin" 0, 



~ ^irfrv^ °> M w) - £& v (0, 0, M w )] - 2n f 2 ^ FV (O, 0, M w ) 
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- nf#2 ( cos 2 6 W + 
n f T 3 



5 sin 6m 



J- 



3 cos 2 

(glAM2 + lg 2 1 AM 1 )+0(g»T 



P{;,. v (0,0,A/ z )- ^e 2 2^ FV (m f ,m f ,0) 



(A.19i) 



24tt wz 3 
3 3 3 

= -|£y^FFs(«7.f,«7,f,«7,l) - ^y^>FFs("7.f,"7,f,"7,2) - -^P F FS ("7-f , 0, m 2 | A.19j) 

1 " , = w(r J1 " J ) + °i» ir, i 



V® = 



2 1 2 



(4tt) 2 



Q-21n2^ +0(g 6 T 4 ) 



(A.19fc) 
(A.19Z) 



^ M = -\g\ cos 2 w X>$ vv (M w , M w , M z ) - ^e 2 ^ vv (M w , M w , 0) 

- -g\ cos 2 6» w [Pllt(M Lw , M Lw , M z ) - P LLT (M W , M w , M z )] 

- \e 2 [V LLT {M Lw , M Lw , 0) - V LLT {M W , M w , 0)] 

- g\ cos 2 ^llt(M Lw , M Lz , M w ) + g\ cos 2 6» W P LLT (M W , M z , M w ) 

- # 2 sin 2 0£> LLT (M Lw , M L7 , M w ) + # 2 sin 2 6> w D LL t(M Lw , 0, M w ) 



^T 3 
16tt 



AM 2 +0(g 5 T 4 ) 



= -2g 2 2 V^ v (M w ) - g 2 2 cos 2 6 W V^ V (M Z ) + ^AA4 2 + 0(/T 4 ) 



-^ 2 ^ V (M W , M w ) - \g\ cos 2 # W 2^ V (M W , M z ) - Je 2 2^ v (M w , 0) 



1 

4 

2 2 2 T 3 



2n 2 T 2 

AA4 2 + ^M W AM 2 + 0(g 5 T 4 ) 



32tt " " ' (4tt) 2 
A. 5. Results for resummed graphs in terms of unresummed graphs 



(A.19m) 
(A.19n) 

(A.19o) 



Pssv(mi, ma, M) =Z>£ sv (mi, m 2 , M) + — (M L - M)T 3 + 0(// 2 T 4 ) (A.20) 

Psvv(m, M l5 M 2 ) = P^w(m, Mi, M 2 ) + 4if 3 (m, M L1 , M L2 ) - 4H 3 (m, M 1; M 2 ) + 0(^T 2 ) 

(A.21) 

P S v(m, M) = vf v (m, M) + J-(M L - M)T 3 + i-(m - m)T 3 + 0(// 2 T 4 ) (A.22) 
2W™ f2 , m n , M) = X^ FV (m f2 , m fl , M) + -j-(M L - M)T 3 + 0(g 3 T 4 ) (A.23) 

Z47T 
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^FF S (m f2 ,m f i,m) =^ FS (m f2 ,m f i,m) + ^(m-m)T 3 + 0(// 2 T 4 ) (A.24) 

Pvw(Mi, M 2 , M 3 ) = V$ vy {M u M 2 , M 3 ) + ^(M L i + M L2 + M L3 - Mi - M 2 - M 3 )T 3 
+ ^ X] [P LL t(M l1 ,M L2 ,M L 3)-Pllt(M 1 ,M 2 ,M 3 )]+0(^ 3 T 4 ) 

perms. 

(A.25) 

The sum above is over all six permutations of (Mi, M 2 , M 3 ). 

(Mi + M 2 )(Mi-M 2 ) J 



Pllt(Mi, M 2 ,M 3 ) 



(4tt; 



2 



(Mi + M 2 )M 3 - MiM 2 + 



+ (Mf - 2Mi 2 - 2M|)# 3 (M 1 , M 2 , M 3 ) 



(M 2 — M 2 ") 2 

+ 1 ' m2 2J [H 3 (M 1 , M 2 , M 3 ) - if 3 (Mi, M 2 , 0)] 

(A.26) 

P w v(M) = P^ V (M) - -L(M L - M)T 3 + 0(^ 3 T 4 ) (A.27) 



V vy (M u M 2 ) =P$ v (Mi, M 2 ) + -^(M L i + M L2 - Mi - M 2 )T 3 

1D7T 



4T 2 

[(M L i - Mi)M 2 + Mi(M L2 - M 2 )] + 0(s 3 T 4 ) 



(A.28) 



Exact results for unresummed graphs in terms of I, H and L 

This section is included only for completeness. The expansions for the V s that give 
our final result for the effective potential through 0(g 4 T 4 ) are given in the next section. 
We haven't bothered to keep track of the Hff terms in the fermion contributions below 
because these do not contribute to the potential at 0(g 4 T 4 ). In the case of Dffs, these 
terms in fact depend on whether the S is a scalar or a pseudo-scalar — a distinction we are 
otherwise able to ignore. 

V% sv ( mi ,m 2 ,M) = I(M)[I( mi )+I(m 2 )} - I( mi )I(m 2 ) 



+ (M 2 - 2m\ - 2m 2 2 )H(m u m 2 , M) 

+ K M2 ml) l^M) - 1(0)] [/(mi) - I(m 2 )} 

+ ^^PK, ma, M) - H( mi ,m 2 , 0)] 
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(A.29) 



V^ vy (m,M 1 ,M 2 ) = {^[I(Mi) -/(0)][J(M 2 ) - J(m) 

9 

m 



^^[/(MO - J(0)][J(M 2 ) - 7(0)] + (5 - 4e)tf (m, Mi, M 2 ) 



+ 
+ 



(A.30) 



2M 1 2 Mf 
(M 1 2 -2m 2 ) r - 



M| 



[if (m, Mi, M 2 ) - 2/f (m, Mi, 0) + H(m, 0, 0)] 

[H(m,M 1 ,M 2 )-H(m,M 1 ,0)}\ + {Mi <-> M 2 } 



P* v (m, M) = -2(3 - 2e)/(m)/(M) 



(A.31) 



2)£ FV (m f2 , ™fi,M) = 2(1 - e)[/ f (m f2 ) + 7 f (m f i)]/(M) - 2(1 - e)/ f (m f2 )/ f (m {1 ) 

,2 * \ (A-32) 



(m| -rofij 
M 2 



[if (m f2 ) - 7 f (m f i)][/(M) - 7(0)] + 0(g 2 T 2 H s ) 



£^ FS (m f2 , mfi, m) = 2[/ f (m f2 ) + J f (m fl )]/(m) - 2/ f (m f2 )/ f (m fl ) + 0(g 2 T 2 H ff ) (A.33) 



^ vv (Mi,M 2 ,M 3 ) = i ^ {l^*' ' )- l^Sl^ )] 2 -^'^) 

perms. 2 3 



4 M 2 2 Mg 



+ 



+ 



+ 



+ 



Mf 



(-4 + 4e)^ + (--4e 



M 2 M 2 1 Mf 



M| 



2 M 2 M 2 



7 „ \ ^9 . \ (Ml Mf 



2 2e J + I 2 4e J I M| M 2 J 4 M 2 M| 



H(M u M 2 ,0) 
I(M 1 )I(M 2 ) 



- -4 

2 e y vm 2 m? 



M 2 M 2 ^ 1 M 2 4 



Mf 



2 M 2 M 2 



9 



/(Mi)/(0) 



(_ 8 + 6e )Mi ! + (4 - 46) ^ + ^ -- + 4e 



x H{M 1 ,M 2 ,M 3 ) 



M 2 M$ 1 Mf 



+ T" 



M| 4 MfMf 



(A.34) 



The sum above is over all six permutations of (Mi, M 2 , M 3 ). 



V vvV (M) = -l/(M)/(0) - ^M 2 /f(M, 0,0) 



(A.35) 



£>£ v (Mi, M 2 ) = 2L(Mi, M 2 ) + (-14 + 20e)I{M 1 )I{M 2 ) (A.36) 
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A. 7. Expansion of results for unresummed graphs 



4 , MT 3 M 2 T 2 

V f ssv (m u m 2 , M) = — — + 

+ 0(/ /2 T 4 ). 



11 1 



2^ vv (m,M l5 M 2 ) = 



T 5 



12 e 12 



5 1 



+ ( ^ ) +ln ( ) + ^ + ^c b -4ch 



J"2 



3T 



M 



1 1 



(4tt; 



2 e 



A 2 



^2 



M 



+ + 10 + 



Mi 
Mf 
M 2 



M| 
M| 



In 



3T 

m7 



M| 
+ M 2 



M 2 T 2 



X>^ FV (mf 2 ,mfi,M) = 



11 1 
h - In 

2 e 2 



^2 



+ 



+ 



M 2 T 



(4tt) 2 
11 1 



2t^2 



6 e 6 



1 1 

'11 



A 2 



J In 

'3T 

,M2 

MT 3 

2 



2 
3T 



Mo 



Mi +M 2 
+ 0{g 1 ' 2 T 2 ) 

mT 3 



(A.37) 



(A.38) 



+ 



In ( ^ ) + 4 + ^c F + ^ln2 



^2 



8tt 
1 1 



+ 0(// 2 T 4 ) (A.39) 



1 



1 1 



-- + - In ^ c B - - In 2 



^ / \ (m 2 +m 2 )T 2 

^FFs( m f2> m fl> m ) = U~Y2 



(4tt) 2 
+ 0(// 2 T 4 ) 



T 2 7 6 3' 
11 1 



2 6 
MT 3 



+ 



24tt 



+ 0(g 3 T*) 



4e 



In 



P$ VV (M,M,M 3 ) 



^2 r 



(4tt)' 



(2M 2 + M|) ( - 



Ji2 



61 1 61 



1 1, „ 
+ -c F + -ln2 



+ 



(A.40) 
mT 3 



24n 



(A.41) 



-In 



H MM 3 Mo H i 

6 3 3 2M 



24 e 24 

M 3 



+ 



13 13 



T 2 8 



f 12 Cb + 2Ch 



M 3 4 
2M S 4M 2 



+ -4M 2 + 9M; 



M 4 
2Mf 



M| 



/ o o 4M 4 Mo 6 \ / 3T \ 

+ (-12M* - 17Mf + + ^ J in (^-^ j 



4M 3 4 

m 2 ~~ 2M 4 y 

Ml \ 



In 



3T \ 
M + M 3 ) 



+ 



M 4 , /3T\ M| , ,. 

In — H ^ In . 

M y 4M 4 V M 3 7 
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A. 8. Some useful limits 
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Appendix B. Derivation of Some High- Temperature Expansions 

Expansion of 

In this appendix, we shall derive the leading 0(T 2 ) term in the high-temperature 
expansion of Hq {jxifz > wifi ; ttl) , defined by ( |5.12|) . Because two of the three Euclidean 
momenta P, Q, P + Q are fermionic, -f/g- does not diverge in the infra-red if we set all the 
masses to zero. So the leading 0{T 2 ) term is simply 

r(i) /-(f) i 

if ff (m f2 ,m fl ;m)=^ ^ ________ + O(mT). (B.l) 

The result found in this appendix is that 

/■(f) 1 

t t p^PTW = m (B ' 2) 

in dimensional regular ization and so vanishes when e — > 0. 
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We begin by rewriting the integral in (|B.2|) as 



(2^)3-^(2^)3-^ (27r) 3- a « g (P + g + 
oqq oaa even v ' K ' K ' ,„ „. 

n j k (B.3) 

1 

X 



[(httT) 2 + p2] [(jTrT) 2 + g 2 ] [(^ttT) 2 + r 2 ] ' 
Using the symmetry of this expression, one may replace 

odd odd even \ any any any even even even/ 

where we will restrict the triple sum on the right-hand side to exclude the case (n,j, I) = 
(0, 0, 0) where all the frequencies are simultaneously zero. By scaling all three-momenta 
by a factor of 2 in the first term on the right-hand side of ( |B.4| ), we obtain 



EEE-^EEE- cm) 



odd odd even even even even 



So we may relate Hq to something resembling the bosonic function H: 

(f) r({) i r r i - s s 

U e[4\n2 + 0(e)]±± ' 9 ° 2 - (B.6) 



P 5q P 2 Q 2 (P + Q) 2 1 Tp Jq P 2 Q 2 (P + Q) 2 

The result will be non-zero when e — > only if the bosonic integrals on the right-hand 
side give a divergent contribution of order 1/e. To relate this result more directly to the 
bosonic function H(m), we need to temporarily put in an infrared cut-off m: 

1 — S Po S qo _ 



p%p 2 Q 2 (P + Q) 2 

hm \H(m)-T 2 1 d P d Q ' 



m^o[ K ' J (27r) 3 - 2e (27r) 3 " 2e (p 2 + m 2 )(q 2 + m 2 )[(p + q) 2 + m 2 } 

It is straightforward to compute the 1/e term of the above three-dimensional integrals 
using standard techniques for loop integrals. One finds that it exactly cancels the 1/e 
piece of H(m) shown in (|3.17| ), from which (|B.2|) then follows. 
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B.2. Expansion of L{mx 1 m2) 



Working from the definition ( |6.3| ) of L(mi, mi 
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An expansion for the remaining integral may easily be found by the method of ( |3.32| ): 
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The high-temperature expansion of L(mx, m^) is then 
L(mi, 7712) = const. — - — — — — — (mf + m 
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Appendix C. Criticism of the Super-Daisy Approximation 

A previous attempt to compute corrections to the ring-improved one-loop potential 



in ref. |14j has relied on the super-daisy approximation to simplify the calculation. In this 
appendix, we shall explain why the super-daisy approximation does not correctly reproduce 
the leading corrections to the one-loop potential. The super-daisy approximation is an 
approximation to Dyson's equations where the effective masses of particles are derived at 
one loop (to any desired order in M/T) from diagrams like Figs. |l| and [5], and then those 
masses are used self-consistently to improve the propagators used to derive the effective 
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masses in the first place. The result is a set of coupled equations which may be solved for 
the effective masses. In pure scalar theory, for example, the equation is given schematically 
by fig. [2^ (ignoring renormalization and counter-terms): 

m ^ m2 A4 P i4^-\^iw^- < ai > 

Computing the one-loop potential with these masses is then equivalent to some approxi- 



mation to the set of super-daisy diagrams, such as shown in fig. 29 



Super-daisy equations such as (|C.1| ), however, do not give adequate approximations 



to these graphs. (|C.1|) should really read 

n«) » m + \xf p - \^f p [pa + n{mp ; Q)2 + n(p + q)] , (c.2) 

where II (P) is the one- loop self-energy. The approximation made in ( |C.1| ) was to replace 
n(P) by 11(0) inside loops. Thus, fig. has been approximated by first computing the 
smallest (outer-most) loops in the approximation that no momenta flows into them, then 
computing the next smallest loops in the same approximation, then the next, and so forth. 
At each stage, the approximation is valid only if the momentum flowing through a loop is 
small compared to the momenta of all the smaller loops that decorate it. And therein lies 
the rub. Once quadratically divergent pieces have been accounted for (say, by normal ring 
resummation as opposed to super-daisy resummation), then most loops are dominated by 
momenta of order M and there is no hierarchy of momenta in a graph. So, for example, the 
0(g 2 M 2 T 2 ) contributions of fig. 20m cannot be correctly computed in an approximation 
that replaces one loop by 11(0). 

The problem can be even more egregious if care isn't taken to avoid mishandling hard 
(P ~ T) thermal loops. Consider the particular contribution to the super-daisy improved 
one-loop potential shown on the right-hand side of fig. [30[ The super-daisy approximation 
instructs us to first evaluate the small loop assuming Q = 0. This small loop is then 
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dominated by loop momenta P of order M. When we replace the result II (P) of the small 
loop by an effective mass insertion 11(0), the large loop becomes quadratically divergent 
and so is dominated by momenta Q ~ T. But this is exactly the opposite hierarchy of 
momenta from the one needed to justify the super-daisy approximation. Ref. [14|] sidesteps 
this disaster by computing the derivative of the effective potential rather than the effective 



potential itself. In the tadpole graph of fig. gl| the large loop never has momentum T when 



it is decorated by smaller loops. But this does not save the super-daisy approximation. The 
large loop is now dominated by momenta Q ~ M, and the O(MT) corrections from the 
small loop are also dominated by momenta P ~ M; the necessary hierarchy of momenta 
is still absent. 



As discussed in section |4.2| , there are a few diagrams where there is a hierarchy of 
scales when the scalar mass is much smaller than the vector mass. The resummation 
discussed in that section does not benefit from the full super-daisy approximation of (|C.1| ) 
and was in any case irrelevant to our final results for the leading corrections to the one-loop 
potential. 



In section \i.4[ we discussed how the most important terms for shifting the VEV from 
its one-loop value are those that depend on logarithms of masses. The super-daisy ap- 
proximation will never generate such terms because it always reduces multi-loop diagrams 
to a hierarchy of one-loop diagrams (with some number of mass insertions) and the re- 
sult I(M) for a one- loop diagram does not contain any InM terms. [See ( |3.12[ ).] Such 



terms do appear in many of the expressions in ref. [I3J], but this is because there I(M) is 
routinely split into separate zero-temperature and temperature-dependent pieces. Though 
each piece depends on InM, the sum does not. Ref. [fnj does consider one diagram in 
addition to the super-daisy one-loop tadpole, and this one diagram will indeed generate 
some InM terms. However the other InM terms from all those diagrams subsumed by the 
super-daisy one-loop tadpole (such as the derivative of fig. |20|m, for example) are lost. 
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Figure Captions 

Fig. 1. Quadratically divergent loop giving rise to the thermal mass. 

Fig. 2. A generic example of a one-loop ring (daisy) graph. 

Fig. 3. Resummation of the propagator in the ring approximation. 

Fig. 4. Qualitative form of the potential as a function of temperature during a first-order 
phase transition. The potentials have been normalized to be zero at the origin. 

Fig. 5. Dominant contributions to the scalar thermal mass in gauge theories with 
fermions, in addition to fig. |l|. 

Fig. 6. Generic one-loop ring diagram in a gauge theory with fermions. 

Fig. 7. 2-loop contributions to the potential of the scalar toy model. 

Fig. 8. 2-loop contribution involving the thermal counter-term. 

Fig. 9. Cancellations of resummation and thermal counter-terms. 

Fig. 10. Resummation of fig. [7]b. Double lines are resummed, "heavy" denotes ko ^ 
modes and "zero" denotes k = modes. 

Fig. 11. Relation between a heavy contribution and the unresummed graph. 

Fig. 12. Rewriting of fig. [TIJ 

Fig. 13. Resummation of fig. 0a. 

Fig. 14. Mixing between the polarization of and the unphysical Higgs boson in a 
background field (f>. 

Fig. 15. 2-loop diagrams in the Abelian Higgs model. 

Fig. 16. A generic multi-loop diagram. 

Fig. 17. Generic once-iterated ring diagram for the Higgs- loop contribution to the effective 
potential. 



Fig. 18. Resummation of fig. |15|a. 



Fig. 19. Two-loop diagrams involving fermions. 

Fig. 20. Additional graphs needed for a non-abelian theory. 



Fig. 21. Resummation of fig. |20|m. 



Fig. 22. Rewriting of fig. |2l]. The L and T lines represent the longitudinal (Aq) and 
transverse polarizations of only the po — contributions. 
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Fig. 23. 2-loop graphs for the Minimal Standard Model, "q" indicates all quarks other 
than the top quark. 

Fig. 24. Relative shift A</>/(/> of the asymmetric vacuum at the critical temperature plotted 
vs. the Higgs mass. The solid line comes from using the full formula for the 2- 
loop potential derived in this paper; the dashed line should approximately match 
it up to yet high-order corrections in the loop expansion. The high-temperature 
expansion used in this paper is valid provided the dotted line (described in the 
text) is small compared to one. 

Fig. 25. The ratio 4>/T at the phase transition. The dashed line is the one-loop result and 
the solid line is our two-loop result in Landau gauge. 



Fig. 26. Same as fig. || but for m t =180 GeV. 

Fig. 27. Same as fig. [25| but for m t =180 GeV. 

Fig. 28. The one-loop approximation to Dyson's equation for the self-energy. 

Fig. 29. A generic super-daisy diagram in the pure scalar theory. 

Fig. 30. A disastrous term in the super-daisy resummation of the one-loop vector diagram. 

Fig. 31. The tadpole version of fig. |3(], which is still problematical. 
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